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An Integrated Approach to the Determina- 
tion and Selection of Lunar Trajectories 


Walter C. Nelson* 


\bstract 


An analytic approach to the determination and selection of 
allistic lunar trajectories based on the restricted three-body 
roblem is presented. A technique for relating a geocentric 
llipse to a selenocentric hyperbola is shown. The technique 
sads to inter-relationships of major design areas including 
widance and control, propulsion, tracking and communica- 
ions, and logistics. Particular emphasis is placed on relating 
ehicle energy to selenographic coordinates for lunar landing 
nd takeoff or for entering and leaving lunar orbits from which 
he “‘corridor-ring” concept is developed. In addition an analy- 
is is conducted that defines the optimum point for entering 
nd leaving lunar orbits or for lunar landing and takeoff. 


I 


ntroduction 


An approach to ballistic lunar trajectories using 
<eplerian orbits as a basic ingredient is presented. A 
eocentric ellipse with earth radius as perigee is matched 
ometrically and dynamically to the asymptote of a 
lenocentric hyperbola. For a lunar landing this hyper- 
la intersects the lunar surface at the landing site 
d for transferring to a lunar orbit, the point of trans- 
ris at periselenum. A proper match of these Keplerian 
its is accomplished with the aid of the energy equa- 
on. The energy equation can be applied to translunar 
hicles to determine the periselenum velocity decre- 
ent required for landing or orbiting. By a reverse 
ocess the velocity increment needed for transearth 
ght can be determined. The fact that the moon 
oves about 3500 ft/sec relative to earth leads to 
sitive energy selenocentric vehicle orbits from which 
try and exit corridors are derived. 

Corridors are defined by the volume swept out when 
tating a segment of a hyperbola about an axis paral- 
I to the asymptote and passing through the moon’s 
nter. The locus of the periselenum point for a given 
sitive energy lunar orbit determines rings for entry 
or exit from a closed lunar orbit when periselenum 
greater than the moon radius. Lunar orbits with the 
me positive energy are used to define rings for lunar 
nding and takeoff when periselenum is less than a 
oon radius. Rings are selenographically fixed if the 
all angular motion of the moon relative to the earth- 
oon line is neglected. The direction of the corridor- 
g axis which becomes fixed when time of flight and 
ajectory inclination are specified defines the avenues 
a vehicle approaching or departing the moon. 
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In addition to the corridor-ring concept, an analysis 
is conducted showing that the optimum point along a 
hyperbolic orbit for entering or leaving a circular lunar 
orbit is at periselenum. 

The ring-corridor concept is an aid to the lunar 
trajectory problem that provides some simple tools 
leading to an appropriate selection of ballistic trajec- 
tories for a given mission and in addition provides an 
integrated outlook to problems relating guidance and 
control, propulsion, communication and_ tracking, 
and logistics. 


Analytical Model of the Earth Moon System 


An approach to lunar trajectories using Keplerian 
orbits as a basic ingredient is presented. A geocentric 
ellipse with earth radius as perigee is matched geo- 
metrically and dynamically to the asymptote of a seleno- 
centric hyperbola. The matching of the ellipse and 
hyperbola provides an analytical technique for ob- 
taining initial conditions for ballistic translunar and 
transearth flight. With the initial geocentric or seleno- 
centric velocity and the corresponding geographic or 
selenographic coordinates, a high speed computer can 
then be used to “zero in”? on a precision trajectory. 
The precision trajectory is the result of integrating the 
vehicle motion in some appropriate mathematical 
model of the earth-moon system. 

The analytical model of the earth-moon system used 
for the precision trajectory calculations consists of 
the earth and moon revolving in elliptic orbits about 
the barycenter with the barycenter in turn moving in 
an elliptic orbit about the sun. The moon is assumed 
to keep the same side facing the earth at all times. 

The moon’s orbit can be described relative to a direc- 
tionally fixed coordinate system with origin at the cen- 
ter of the earth by two geometrical parameters 


eccentricity em = 0.0549005 
semi latus rectum p,, = 1.2585193 X 10. 


The mass characteristics of the earth and moon corre- 
sponding to the eccentricity and semi latus rectum can 
be represented by the respective products of their 
masses and Newton’s gravitational constant 


ue = 1.40769 x 10° ft.’/sec” 


1.7299 x 10” ft.*/sec”. 


K m 


The Journal of the Astronautical Sciences 33 


CORRIDOR 


4 


AA 


AXIS OF 


Retgae | REVOLUTION 


MOON 7 


ASYMPTOTE 


PERISELENUM 


FOCUS 


RING 


Fie. 1. Segment of hyperbola used to generate corridors 
and rings. 


The radii of the earth and moon are assumed to be: 
re = 20.92 X 10° ft. 
tm = 5.7024 X 10° ft., 


respectively. 

The constants used to represent the moon’s orbit 
and the mass characteristics of the earth and moon 
are regarded as being typical but not necessarily con- 
sistent with the actual values. 

The foregoing mathematical model is used to gener- 
ate ballistic trajectories for translunar and transearth 
flight. The individual mass characteristics of the earth 
and moon and the geometrical arrangement of the 
earth and moon are in turn used in Keplerian orbit 
mechanics to estimate initial conditions for the trajec- 
tory integration. 


Six Steps for Orbiting the Moon 


A trip to the moon incorporating lunar orbiting can 
be visualized as taking place in six steps. Assume a 
trajectory coplanar with the lunar orbit plane. 

1. Injection on to a geocentric ellipse 

2. Transfer to a selenocentric hyperbola—‘‘corridor”’ 

3. Enter lunar orbit or land on the moon—‘ring”’ 

4. Takeoff from the moon or leave lunar orbit— 

“ring”? 
5. Injection on to a selenocentric hyperbola- 
“corridor” 

6. Transfer to a geocentric ellipse 
The primary idea contained in these steps is in the 
corridor-ring concept. The surface of the corridor is 
formed by rotating a segment of a hyperbola about an 
axis passing through the center of the moon, one 
asymptote being parallel to the axis of rotation. The 
center of the moon is assumed to be at the focus. The 
locus of the periselenum point on the hyperbola forms 
the ring on the surface of the moon or on the surface 
of a sphere concentric to the lunar sphere. The method 
for generating the corridor and ring is illustrated in 
Figure | 


The geocentric ellipse of step 1 has a perigee near 
the surface of the earth with apogee extending to the 


: ae : : Spas 
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Fig. 2. Illustration of the six steps for orbiting or landi 
on the moon and returning to earth. 


near side of the moon or beyond the moon. A vehicle 
traveling to the moon moves away from the earth 
Figure 2—on this eu When it has reached 
vicinity of the moon’s orbit the vehicle velocity = | 
be low in relation to the moon’s velocity. This i 
effect coupled with the gravitational attraction of th 
moon, transfers the vehicle on to a selenocentric hyper 
bola placing the vehicle on the surface of a corrid 
As the moon moves toward a collision with the vehie 
the vehicle is in effect moving along the surface of th 
corridor toward periselenum. At periselenum th 
vehicle is on the entry ring and application of an appro 
priate impulse puts it in a lunar orbit. Orbiting con 
tinues until return to earth is_desired. 

In general the point for entry of the lunar orbit wil 
not be coincident with point for exit. However, one 
during each circumnavigation the vehicle passes over 
point on an exit ring where it is suitable for initiatin 
escape to return to earth. Then application of the ap 
propriate impulse causes the vehicle to move away 
from the moon on a selenocentric hyperbola. Applica 
tion of the escape impulse can be regarded as permit 
ting the moon to move away from the vehicle whie 
tends to maintain a fixed position relative to the eart 
As the distance between the moon and vehicle int 
creases, the vehicle velocity decreases and it transfer! 
on to a geocentric ellipse and descends to earth. Periges 
is assumed to be at the earth surface. 

A quantitative description of the foregoing luna! 
mission can be formulated using the mechanics 0} 
Keplerian orbits. To do this a backward approach will 
be used and it will be assumed that the vehicle is in ¢ 
lunar orbit and it is desired to bring it to earth ballis: 
tically. Two assumptions will be made concerning th 
lunar orbit for purposes of illustration, but these as 
sumptions need not apply generally. First assume thai 
the orbit is circular at zero altitude. This means thi 
orbit velocity is given by (tm/?m) '? Next the line 
nodes formed by the intersection of orbit plane and thi 
lunar orbit plane is perpendicular to the earth-moo1 
line. This line is coincident with the velocity vector ¢ 
the moon when the moon is at apogee and perigee. 

Ballistic return to earth cannot be accomplished by 
providing the vehicle with a velocity only sufficient t 
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Fic. 3. Illustration of relation between geocentric ellipse 
1 corridor for a minimum energy mission. 


sape the moon (2um/rm)”. This can be seen ex- 


ul. BK 
Sac Sie (1) 


scape speed implies zero total energy. Therefore 
hen the vehicle has escaped the moon and the dis- 
nee from the moon is large, the relative velocity 


directionally fixed earth centered coordinate system 
is about 3500 ft/sec and essentially perpendicular to 
3 earth-moon center line, it follows that the vehicle 
locity relative to earth following lunar escape must 
be roughly the same as the moon’s. That is, the veloc- 
4 of the moon is inherent to the velocity of the vehicle. 
Consequently a geocentric vehicle orbit with dimen- 
jons similar to the dimensions of the moon’s orbit 
will result. This problem can only be avoided by es- 
caping the moon at hyper escape velocities. 
The increment of velocity above escape needed to 
bring a vehicle to earth can be deduced from equation 
(1) by rewriting it in the form: 


(ve = Av)’ aad S oe 5 Veo (2) 


where Av is the excess vehicle velocity above the es- 
eape velocity v, and v.~ is the velocity relative to the 
moon for large values of r. It follows that 

Av = (v¢ + ve)? — % - (3) 
This equation states that if the velocity at r — © is 
cnown, the amount Av needed to be applied at the 
unar surface to produce v.. can be determined. 

A low energy transfer from the moon to earth can be 
ised to illustrate the application of equation (3). 
Suppose it is desired to transfer a vehicle from a 
elenocentric escape orbit on to the apogee of a geo- 
entric ellipse. The apogee velocity relative to earth is 


tiven by: . 
D) 1 1/2 
VG = E (; a ‘) ) (4) 


vhere r is the distance from the center of the earth to 
he surface of the corridor and 


Co 4(r =P (ao (t y 


Cr 
—S 


where r, is the earth radius—see Figure 3. It is desired 
that v.. be equal to the difference between the apogee 
velocity of the geocentric ellipse and the velocity of the 
moon. 


Vn = Um — Va (6) 


where v, is the magnitude of the moon’s velocity. 


The relation of v. to the geometrical and dynamical 
properties of a selenocentric hyperbola with periselenum 
at the lunar surface can now be derived. From the 
Keplerian orbit equation 


Di Pm(1 ie é), €7') 


where p is the semi latus rectum of the hyperbola, e is 
the eccentricity, and r,, is the moon radius. Equations 
(3) through (7) can be used to derive expressions for 
the geometrical description of the escape hyperbola in 
terms of the required escape speed. 


2 

a —1 (8) 
2 

pe eae (9) 


Having obtained the geometrical properties of the es- 
cape orbit the dimensions of the corridor and ring fol- 
low and these in turn are related to the velocity applied 
for escape. The corridor radius is as follows: 


Pm 


Di Oe 
E i (5 + =) | 


The quantity 6 is called the guidance parameter ac- 
cording to Mickelwait in reference 3. Relations be- 
tween Av, e, and p are given in Figure 4. The ring size 
is specified by the cone angle subtended at the center 
of the moon: 


b= (10) 


= 180° — cos — (-!). (11) 


1.0 2.0 3.0 4.0 5.0 6.0 7.0 


e- ECCENTRICITY as 
p-SEM! LATUS RECTUM~% OF r,, (IO °) 


Ov-EXCESS VELOCITY IN % OF PERICENTER ESCAPE SPEED 


Fig. 4. Semi latus rectum and eccentricity of a hyperbola 
vs. excess velocity applied at pericenter. 


The Journal of the Astronautical Sciences 35 


TABLE 1 


Comparison of the low energy mission obtained from 
Keplerian orbit calculations to integrated 
trajectory calculations 


Parameter Keplerian Integrated Trajectory 
Corridor radius 2.85 2.93 
(moon radii) 
Ring size—p 39 .5° Bien 
Av ft/sec 526 499 
Geocentric ellipse 966 967 
eccentricity 
Time of flight 3.098 Apogee to | 4.531 Periselenum 
(days) perigee to perigee 


Equation 10 shows that as the injection speed v, + 
Av increases the corridor radius gets smaller and equa- 
tion 11 indicates that the ring size increases because e 
increases with escape speed. Some limitations on Av 
will be discussed in subsequent paragraphs. 

Some specific results for the low energy mission have 
been tabulated in Table 1. The values tabulated under 
the integrated trajectory were obtained using the 
initial conditions from the trajectory computation in 
conjunction with equations (10) and (11), except for 
the geocentric ellipse eccentricity and time of flight. 
These were extracted directly from the computed tra- 
jectory. It is seen that the corridor radius is small in 
comparison to the distance to the moon. Because of 
this it can be shown that the Av needed to orbit the 
moon or escape a lunar orbit on the corridor surface 
farthest from the earth will differ a negligible amount 
from the given value. This is because the apogee 
velocity of a geocentric elliptic orbit changes only a 
sight amount for changes of apogee over the dimen- 
sions of the corridor diameter, which in effect maintains 
Vv. constant. On the basis of this result the rings and 
corridors are defined to be circular. It is of interest that 
subsequent trajectory computations showed that get- 
ting off the moon on a trajectory that moves across the 
back side of the moon, over the top, underneath, and 
across the visible face took virtually the same amount 
of energy. These results can be of interest to communi- 
cations and logistics of a mission design in addition to 
the propulsion and guidance. 

If it is assumed that the moon is in a circular orbit, 
the trajectory of a vehicle traveling to the moon can 
have symmetry with a trajectory that goes from the 
moon to earth, as illustrated in Figure 5. A round trip 
trajectory is plotted showing relative positions of the 
moon and vehicle as functions of time for a low energy 
mission. The coordinates are earth centered and direc- 
tionally fixed. Instantaneous transfer from the entry 
ring to the exit ring is assumed along with application 
of the appropriate velocity increment for entry and exit 
from the lunar orbit. The inset shows the vehicle tra- 
jectory in the vicinity of the moon in a directionally 
fixed earth centered coordinate system. 
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Fia. 5. Low energy round trip to the moon. Trajectory co- 
planar to moon orbit plane. 


Lunar Landing and Takeoff 


Lunar trajectories applicable to missions incorporat- 
ing landing and takeoff can be developed with the aid 
of the corridor-ring concept. It was indicated that the 
energy relative to the moon required for a moon t¢ 
earth transfer on the surface of a given corridor is 
essentially constant. This was because the surface of the 
corridor represented the locus of asymptotes formed by 
rotating a given hyperbola about an axis of revolutior 
through the center of the moon. Consider any seleno: 
centric hyperbola with this same total energy relative 
to the moon but suppose periselenum is smaller, that is 
under the surface of the moon. This means a smalle! 
diameter ‘‘corridor’”’ can be generated with a hyperbols 
intersecting the lunar surface as shown in Figure 6 
It can be deduced from Equation (3) that constant A 
is associated with any of these hyperbolas at the poin 
of intersection with the lunar surface. This is becausi 
v.». Was assumed to be constant and », is necessarily con 
stant for all points on the lunar sphere. The distinctio1 
between any two given landing or takeoff points o1 
the lunar sphere in relation to a given Av is the angle a 
which the hyperbola intersects the lunar sphere. Thi 
angle can be called the launch or landing angle. Figur 
7 shows how this angle is related to points on the luna 
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_ Fic. 6. Illustration of corridors and rings for lunar landing 
and takeoff. 
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Fic. 7. Lunar landing and takeoff angle vs. selenographic 
longitude in the selenocentric vehicle orbit plane for the low 
energy mission. 


equator for the low energy mission. Points on the lunar 
sphere that have the same landing and launch angle 
form circles concentric to the axis of the approach or 
exit corridor. 

In general an approach or exit corridor is associated 
with a ring whose cone angle is less than 180°. There- 
fore, for a given corridor a majority of the lunar sur- 
face can be regarded as exterior to a given ring. And 
because constant values of v., and Av are assumed to be 
associated with any approach or exit hyperbola in- 
terior to the corridor with asymptote parallel to the 
corridor axis, it follows that most of the lunar surface is 
available for either landing or takeoff using a single 
impulse. 


Time of Flight 


To this point discussion has centered on escaping a 
lunar orbit and transferring to a geocentric ellipse at 


apogee. This is regarded as a minimum energy transfer 
and since essentially the complete half of a geocentric 
ellipse is traversed the flight is regarded as taking maxi- 
mum time. Because shorter flight times are of interest, 
particularly for manned missions involving solar flare 
threats, it will be shown how the corridor-ring concept 
can be broadened for reducing flight time. It happens 
that as time of flight decreases, the corridor axis tilts 
toward the earth moving the ring toward the invisible 
side of the moon. This can be seen quantitatively by 
noting that when leaving the moon the vector d.. should 
have a component directed toward the earth if flight 
time is to be shortened. Therefore equation (6) takes a 
vectorial representation: 


De = Vi v, (12) 


where 0 is a velocity associated with some point below 
apogee on a geocentric ellipse as indicated in Figure 8 
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Fic. 8. Illustration of technique for estimating reduced 
flight times. 
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Fic. 10. Vector diagram illustrating vector relationships 
for entering inclined geocentric orbits on the low energy 
mission. 


Using this approach, estimates for the changes in corri- 
dor and ring sizes were made along with changes in 
direction of the corridor-ring axis for missions with 
flight times less than 4.46 days. A high speed digital 
computer was then used iteratively to obtain curves of 
ring size, corridor size, direction of ring-corridor axis, 
and Av for various flight times, Figure 9. The results 
show that fuel penalty for missions under 3.5 days in- 
creases very rapidly. Also, it can be deduced from the 
curve that a ballistic cireumnavigation of the moon 
that intersects the earth-moon line at one moon radius, 
from the center of the moon takes about 5.2 days from 
perigee to perigee. The ring size for this mission is 
p. = 48°. 


Inclined Orbits 


To reduce time of flight, it was seen that the corri- 
dor-ring axis rotated in the plane of the moon’s orbit. 
Tilting the axis out of the moon’s orbit plane causes 
the plane of the vehicle orbit relative to the earth to be 
inclined to the moon’s orbit plane. This is important 
because it permits selection of trajectories that avoid 
the earth radiation belts. Tilting the corridor-ring axis 
out of the moon’s orbit plane means there will be a 
component of 0, that is perpendicular to the moon- 
orbit plane. For a 4.5 day flight the maximum value of 
this component should be equal in magnitude to v, . It 
follows that 


Vo = [(Um — Va COS 7)” + (vq sin 7 (13) 


A vector diagram for equations 12 and 13 corresponding 
to the maximum time of flight mission is shown in 
Figure 10. The orbit inclination 7, results from the 
combination of the perpendicular component and the 
coplanar component, with v, maintaining constant mag- 
nitude for all inclinations. 
A relationship can be derived between the inclination 
I of the corridor axis and the inclination 7 of a geo- 
centric ellipse relative to the lunar orbit plane for a 
4.5 day transfer. From Figure 10 
hatie hee a fees -. (14) 
Um — Va COS 4 
Tilting the corridor axis about 10° produces an inclina- 
tion of 90°. Inclinations greater than 90° are associated 
with retrograde geocentric orbits. Figure 11 gives re- 
sults for the relation between J and 7 as obtained from 
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Fic. 11. Inclination of vehicle orbit relative to lunar orbs 
plane vs. inclination of ring—corridor axis and require 
velocity increment. 
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trajectory calculations and the comparison with equé 
tion 14 is shown. Figure 11 also shows how the & 
increases with 7. For the retrograde orbit with 7 | 
180°, Av = 1031 ft/sec which is a substantial increas 
over the minimum energy requirement of 499 ft/sec. 


Optimum Point on a Hyperbola for Transfer t 
a Circular Orbit . 


An analysis has been conducted which shows thal 
the optimum point for entering or leaving a circulal 
lunar orbit from a hyperbolic approach trajectory is al 
periselenum. On the basis of the analysis it is con 
cluded that lunar missions incorporating circular luna 
orbits should be designed around a trajectory that ha: 
periselenum at the altitude of the desired orbit. 

Suppose a vehicle is approaching the moon on a hy 
perbolic path as indicated in Figure 12 and transfer to 
circular orbit is desired. This means at some poini 
along the hyperbolic path the total energy must be de 
creased to a negative value. The total energy per uni 
mass 1s 


Dy 
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Fie. 12. Hlustration of technique for entering a circule 
lunar orbit from a hyperbolic orbit. 
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Fie. 13. Velocity increment required to enter circular orbit 
‘om a selenocentric parabola vs. distance to moon center and 
ue anomaly. 


The circular orbit velocity is 


othe 


The increment of velocity required to enter a circular 


hit is 
* 
p 


(16) 


a's |] 


bv = [(vsiny)’ + (v cosy — »,)’]}”, (VO 


where y is the vehicle flight path angle at any point 
‘along the hyperbola. Equation (17) can be rewritten: 


, : 2 Fakes 
io = | on + 3H By) (18) 
T fy vie 


Where h is the angular momentum per unit mass. To 
determine if there is some point on the hyperbolic 
path that requires a minimum or maximum 6v to get 
into a circular orbit, differentiate 6v with respect to r 
and set it to zero. It follows that 


2 


fr = fe for ue (6v) = 0. (19) 
bh dr 


Since h’/y is the semi latus rectum of an orbit, equation 
(19) states that for a given hyperbolic orbit the energy 
needed to go into circular orbit is a maximum or mini- 
mum when the vehicle is at the semi latus rectum point 
or when it is at infinity. A curve for the required 
velocity increment versus distance from the center of 
the moon is given in Figure 13 assuming parabolic 
motion with periselenum at the moon surface. It is il- 
lustrated that the velocity increment at the semi latus 
rectum point is a maximum. The required velocity in- 
erement at periselenum is a minimum although it is 
not a mathematical minimum. It is concluded that the 
best point to enter circular orbit from an excess energy 
orbit is at pericenter. 

If the class of all excess energy orbits having the same 
pericenter is considered, a curve of the velocity in- 
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Fig. 14. Velocity increment required to inject into circular 
lunar orbit at pericenter and at r = p. 


crement at pericenter in excess of circular orbit speed 
can be prepared and this is presented in Figure 14. 
Since all the excess energy orbits have the same peri- 
center distance it is convenient to use the eccentricity 
or semi latus rectum as the abcissa. 6v is given in % of 
the circular orbit velocity. For purposes of comparison 
the velocity increment required for entering circular 
orbit at the semi latus rectum point is also presented. 
It is apparent that a substantial fuel saving is to be 
realized by selecting the pericenter point along a hyper- 
bolic orbit for entry departure from a circular orbit. 
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Application of Hill’s Lunar Theory to 
the Motion of Satellites : 


P. Lanzano* 


Abstract 

Hill’s Lunar Theory of Celestial Mechanics is applied to 
the problem of establishing a permanent artificial satellite on 
a periodic orbit around a planet. 

Using a method developed by C. L. Siegel, in his ‘“Vorle- 
sungen tiber Himmelsmechanik,” the Hill’s equations of the 
Lunar Theory are solved to obtain the coordinates of the 
periodic trajectory as Fourier series of the time with respect 
to a rotating system of reference. A recurrent procedure is ob- 
tained for evaluating the coefficients of the series in terms of 
the period of revolution. The Jacobi constant of the motion 
is also expressed as an infinite power series of the period. The 
convergence of such expansions can be ascertained for small 
values of the period. A numerical example for a satellite of 
Venus is furnished. 

An error analysis is undertaken by studying solutions of 
Hill’s equations lying in a neighborhood of a periodic orbit 
and corresponding to the same value of the total energy. The 
coordinates of such neighboring trajectories are determined as 
isoenergetic displacements referred to the intrinsic reference 
formed by the tangent and normal lines at the various points 
of a periodic orbit. This procedure leads to a differential equa- 
tion of the Mathieu type whose solution is obtained as a series 
expansion valid for small values of a parameter. 


Introduction 


The problem of establishing a permanent artificial 
satellite around a planet of our solar system (say Venus 
or Mars) is becoming of actual concern. For many 
reasons, among which are the stability and tracking 
questions, a periodic trajectory is the ideal solution. 

The trajectory that takes the ballistic missile from 
the neighborhood of the Earth (of the order of a couple 
hundred miles above the Harth’s surface) to the vicinity 
of the planet belongs decidedly to the restricted three- 
body problem. 

Extensive numerical studies of such trajectories for 
lunar purposes, whose results are however applicable to 
interplanetary flight, have been undertaken by Russian 
scientists see Ref. [9], but none of these trajectories 
achieves a permanent periodic orbit. It appears, there- 
fore, due to the present limited knowledge, that the 
missile should have a retrorocket which, by slowing it 
down, would ensure its permanent capture by the 
planet. 

However, even within the sphere of action} of the 

* Member of the Technical Staff, Space Technology Lab- 


oratories, Inc., Los Angeles, California. 
t For the definition of sphere of action of a planet see [9]. 
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planet, an approximation of the missile trajectory by 
Keplerian motion does not seem to yield accura 
enough results: various perturbation procedures mu 
be used to obtain, with a greater degree of certaint, 
the position and velocity of the satellite. 

The Lunar Theory of Celestial Mechanics develope 
by G. W. Hill, at the end of the last century, as a pa 
ticular case of the restricted three-body problem ar 
later on elaborated upon by E. W. Brown, appea 
today as the most precise treatment of any perturbatic 
procedure. The latest developments, due to C. 1 
Siegel, even allow the set up of recurrent procedur 
which are highly desirable for numerical calculations t 
means of electronic computers in order to evaluate I 
coefficients appearing in the series representation of tl 
trajectories. The present paper elaborates such pr 
cedures and applies them to the case of periodic orbit 

In the restricted three-body problem, we consid 
two bodies of finite masses moving according to cireul: 
Keplerian motions around their centroid and study tl 
motion of a third body of infinitesimal mass in tl 
plane of the circular motions. The third body is sy 
posed to be attracted by the two finite bodies accordir 
to Newton’s law of gravitation, without affecting the 
relative motion because of its infinitesimal mass. 

In the Lunar Theory, we suppose that one of the ty 
finite bodies has a mass considerably smaller than tk 
other, as might be the case of Venus or Mars wit 
respect to the Sun, and study the planar trajectories « 
a third body of infinitesimal mass which encircle on! 
the body of small mass. 

The periodic orbits of the Lunar Theory can |} 
shown to take into consideration the main solar pe 
turbation on the circular orbits of a two-body proble 
and play the role of the so-called “intermediate 
orbits in perturbation theory. 

In the technical problem of establishing a period 
orbit, three questions belong to the domain of celesti 
mechanics: 

1. The value of the Jacobi constant, in terms of fl 
period of the orbit to be achieved, will allow us— 
ascertain the desired velocity at the time of retrorock 
burnout. 

2. Series expansions for determining the position a1 
velocity of the missile along a periodic trajectory; tl 


S$ necessary for tracking purposes and will also enable 
1S to determine values of the acceleration to which 
the satellite will be subject, information of actual con- 
ern in the case of a manned space station. 
3. Procedures to study trajectories, which can occur 
in the neighborhood of a periodic orbit, in order to 
ealize what will happen in case the periodic traj ectory 
cannot be achieved due to uncertainties in some engine 
parameters. 
The first part of this paper deals with the series 
expansions of periodic trajectories and expression of 
e Jacobi constant in terms of the period. 
_ The assumptions which constitute the Lunar Theory 
reduce the equations of motion for the restricted three- 
ody problem to a particular set of differential equa- 
tions known as Hill’s equations. 
General existence theorems for solutions of differ- 
ential equations assure that under certain conditions 
satisfied by Hill’s equations, the solutions can be 
expressed in terms of infinite power series in the vari- 
bles, within a certain interval of convergence. The 
Becedure consists in assuming, as solutions of Hill’s 
equations, appropriately chosen infinite power series 
two variables with undetermined coefficients; the 
two variables are moreover considered solutions of an 
auxiliary differential system which is known to possess 
periodic solutions. In trying to satisfy Hill’s equations, 
a recurrent procedure arises whereby the coefficients 
of each power can be uniquely determined, once the 
coefficients of lower powers are known. 
_ The absolute convergence of such expansions can be 
established by a “‘majorant’? method which consists in 
comparing the trial series with series having larger 
coefficients and whose convergence can be easily 
ascertained. 

In replacing the two variables by particular periodic 
solutions of the auxiliary system, we obtain periodic 
solutions for Hill’s system with the same value of the 
period. 

By rearranging the various terms of the series, these 
periodic solutions ultimately appear as Fourier series 
of the time; the Fourier coefficients being infinite 
power series of the period, convergent for small values 
of the period. 

The recurrent procedure mentioned before enables us 
(0 determine as many terms as might be required in the 
‘xXpansions of the coefficients of the Fourier series. 
Some of these terms have been obtained by hand com- 
mutation. 

A similar procedure is used to express the Jacobi 
‘onstant of the motion as a convergent infinite series of 
he period. 

A numerical example for a satellite of Venus has been 
vorked out; for this purpose, a program on the IBM 
04 has been used to determine the coefficients of the 
*ourier series up to the twelfth power of the period. 
in idea of such numerical results is furnished by 
Fable 1. 


un 


The last part of the paper deals with trajectories 
lying in the neighborhood of a periodic one and which 
correspond to the same value of the total energy. This 
study can be interpreted either as an error analysis for 
periodic trajectories or might be thought of as furnish- 
ing the secondary perturbation effects due to the Sun. 

Referred to the intrinsic reference constituted by the 
tangent and normal lines at points of a fixed periodic 
orbit, these neighboring trajectories can be specified by 
means of their tangential and normal displacements. 

The normal isoenergetic displacement satisfies, for a 
certain degree of approximation, a differential equation 
of the Mathieu type. 

A modification of a well-known method for solving 
such equations by an infinite cosine series is outlined 
which is very convenient for rapid calculation of the 
coefficients. 

The parameter, which appears in the argument of the 
cosine functions, is also determined by an elementary 
method which, although mathematically rigorous, 
circumvents the question of infinite determinants. The 
paper is meant to be analytical in character and while a 
brief survey of the Lunar Theory must necessarily be 
exposed for reasons of clarity, it emphasizes primarily, 
within the framework of mathematical rigour, new 
numerical methods for determining the essential 
parameters of motion. 


Periodic Orbits 


We shall denote the planet by P; and the Sun by P2 , 
with masses m, and m, respectively. These two celestial 
bodies will be supposed spheres, homogeneous in con- 
centric layers. Neglecting the eccentricity of the 
planet’s orbit around the Sun (in the case of Venus this 
eccentricity is approximately 0.00681, the smallest in 
the solar system), we shall assume P, as describing a 
circular, uniform path around P, with sidereal period 
7*. The missile will be denoted by P;, with zero mass 
(m; = 0) and will be supposed to move in the plane of 
the circular motion, subject to the Newtonian attrac- 
tions of the planet and Sun without, however, dis- 
turbing the Keplerian motion of these two bodies. The 
study of the motion of P; (infinitesimal body) con- 
stitutes the so-called restricted three-body problem. 

The following fundamental units are chosen: 


Ry = distance (P:P2) as unit of length, 


m, + ms as unit of mass, 
* . 
as unit of time. 


20 
It follows then, that 


aol ste , with 0 < w < 1 will represent the 
my, + Me 


mass of P;, 


Meo . 5 
1-—p.= “will represent the mass of P2, 


m + Mm 


i ae et 
w* = — will be the unit of angular velocity and the 
Te 
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universal constant of gravitation 
(w*) "Ris 
my + me 
We consider, next, a Cartesian orthogonal frame in the 
plane of the circular motion with origin at the planet 
P,(0; 0) and rotating with constant angular velocity 
equal to one. 

With respect to such reference, the Sun will occupy 
the fixed position P:(—1; 0) and the differential equa- 
tions for the motion 3 P;(x; y) are given by (see 
Ref. [8] page 347, Ref. [6] page 291-95): 


Wgtet+td—-»vt+F, 


will have the value one. 


G= 


a 
p= ty +h, (1) 
F(e;y) = oii! 4 = 


[(a + by + y?| we (a + app ye? 


where dots denote derivatives with respect to the time, 
and letters used as subscripts, partial derivatives with 
respect to the same variables. 

Introducing the complex, conjugate variables (uw; v) 
by means of the regular transformation of coordinates 


1 : 
teagan iy) 


(2) 
eS oy eer eee cai, 
bw 
the system (1) will be taken into 
&é@ = —2iu + A, 
d= 20+ HA, 
eee 2 tei 
A(u;v) = w+ NE + AD (u + ») (3) 
2(1 — p) 


_ 


pee ao Up'!?) 1/2] +- ies 


If the mass » of P; is so small compared to the mass 
1 — yw of P, that we can neglect the positive powers of 
uw’ in the expansion of H( u;v), then we get the Lunar 
Theory provided we limit ourselves to the study of 
trajectories surrounding only the planet P,. This 
assumption is certainly plausible in the case of Mars 
and Venus. For Venus » = 1/408,001 approximately. 
We have then the following: 


28 3 2 
H(u;v) =- = Fu Oo 
a T Gyn 
4 ome Selig Mera ot : si 
plus positive powers of »'’, hence we get Hill’s differ- 


ential equations: 


| oo 


@= —2i4 +-(ut+) -— — 


aa a 
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SS 
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(4) 
b= 26 +5 (ut) = 


‘ay 


We shall seek solutions of Hill’s equations in the form 
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of infinite power series in two variables with - undete 


mined coefficients. iy 
For this purpose let U(t) and V(t) be period: 
solutions of the auxiliary differential system 


U= a; V = —-aV 
t 2 
= + —;3 ~ = —1) 
a Se (UV) 2 ( ? 
and consider the expressions: 
Vacs [psitsey3iae : 
G33 a 
WwW. , = [prey 
Diet ae q 
Here j is any positive integer and k any integer (p 
tive, zero or negative) satisfying the inequality 
= 4 |kl, 
thus, every exponent in (6) shall be positive ani 
moreover 


t 
‘ 
i 


~~ 
iG 


W aaW x= Wosepe 
The solutions we propose to consider are: 
u= UA + die Ain i) 
v= Vil + Die Ag poe 


the Aj.,.’s being the coefficients to be determined ai 
where the summations shall be considered as doub) 
summations, the index 7 running from one to + © a 
for each 7 we must consider the sum of the finite numbe 
of terms corresponding to k such that 4 | k | S$ 37; e 


Be A j-nW 5% = Ay,0W1:0 —- (Ao, aWo,4 
+ Ao,oWe. + AoiWo1) > 


By formal manipulation of the series (7), we shall fin 
the conditions to be satisfied by their coefficients | 
that the (7) be solutions of the equations (4). At 
later stage, we shall discuss the absolute conve al 
such expressions thus justifying the formal operations: 
For sake of conciseness, we introduce the notation: | 


Ae ate hae 
B= — Don AinWir (& 
C= (1— AO = By ee 


By formal expansions of (1 — A)“? and (1 — By 
into power series of A and B, using the first two e 
pressions (8) and (6’), C can ultimately be expresse? 
as a double summation in the W ;.;’s. Ditferentiatin 
(7) termwise with respect to the time and climinatin 
U, V through (5), we get: 


| 


= 4aU"[L + Qi i(2k + 1)AjaW jal | 


— —4aV"[1 + os (2k +- 1L)Aj.nW 5.x] 

: é 7 
ti = (4a)’U[l + Dog (2k + 1) AW 
6 = (4a)°V"[1 + Doge (Qk + 1)?Aj.W;al. 


Multiplying the two Hill’s equations by U’V° a 
U"V" respectively, replacing in them the (7) and (7 


for the unknown functions and their derivatives, and 
aking use of the notations (8); we find that the fol- 
owing expression: 


oan (2k + 1)"AseWsu] — $4 — 3B 
= CF Wi + dine (26 + 1)ApeWysaul (9) 
— #Wo(1 — A) — $W2,101 — B), 


‘must be an identity in the W;.;’s. 

If we denote by Cj, the coefficient of Wj., in the 
expression appearing in the right-hand side of (9), it 
is possible to realize that C;,, turns out to be a poly- 
nomial (i.e., a finite expression) in the A,.,’s with 
mm < 7}. 

This property allows us to write from (9) the recur- 
‘rent linear system: 


[e+ 1 + 2k) "Ag + GAR = Cir 
RO) Age ia (lL = 2k) A = Cpa, 
Ajo = Cj) 
(10) 
whose solution is easily found to be: 

, - alee are 

Ss 4 (4 — 1) ~* 

2 ae = Tg ie eae oe = 

, Ajo = aCj.0. 


It is possible to ascertain that the series in question, 
whose coefficients satisfy the linear system (10) ( being 
; thus formal solutions of Hill’s equations), are abso- 
_lutely convergent for small values of | U | and | V |. 

The proof is based on a ‘“majorant’”’ procedure and is 

omitted from this paper; an idea of such proof can be 
_ obtained from [5] pages 89-92. 

Among the solutions of the auxiliary system (5), we 
shall choose the periodic one 


U =o” exp(at); V = o'” exp(—at) 


1 " (5’) 
4a = Bai = —1), 
oO 
where o” is a real quantity. 
Thus, if we introduce the notation: 
Ap = 14+ D0; Ajo” 
1 
Ay = 20; Ajuo’” (12) 


ties 0, 37-24. | 4.) 


Hill’s equations will have the following exponential 
solutions: 


= 2k de 
ee ye A, exp (+ = it) 


In terms of the Cartesian coordinates in the rotating 
system 


I : 
L= LC + u); ap 5H —u), (14) 


the solutions for the motion of the third body appear 
as infinite series of the Fourier type: 


i ope Dey Aj, cos ( zi : ) 


o 


and (15) 
_ . (2k+1 
y=at op da A; sin ( ms ). 


These formulae represent two families of periodic solu- 
tions for the motion of P; around P,; because the period 
7 = 2mo° is an arbitrary parameter, within the region 
of convergence of the coefficients, and the double sign 
appearing in the second equation corresponds to the 
possibility of a direct or retrograde motion of P3 
around P; as compared to the circular motion of P; 
about P»2 . 

It is easy to realize that these periodic solutions are 
also symmetric with respect to both axes of the rotating 
frame so that only one-fourth of the trajectory need be 
studied. 

The coefficient A, of these series solutions, as defined 
by (12), are power series in the parameter o, con- 
vergent for small values of the parameter and starting 
with terms in o’/*’ at least. By means of the recurrent 
procedure, illustrated above, it is possible to find as 
many terms as might be necessary in the expression 
for the A;’s even by hand computation. 

The following approximations have been found by 
hand computation up to the fourth power of o”: 


20 fier '6 4 9 19.565 12 
Ao = 1 — 30 + ree sro + e33087 
3 6 ie) 31 12 
Ay = 167 +- go + 967 (16) 
@) 
fe Tog 157,12 
ASI 62. 8 2887 
Dye 
A» sm 560 ) 


all other coefficients start with higher powers of the 
parameter. 


The Jacobi Constant in Terms of the Period 
Hill’s differential equations admit the first integral 


eer 2 2 
“a —=(+ 0) —- ae > —J, (17) 
4 (wv) 1/2 
where J is the so-called Jacobi constant. 
In terms of the Cartesian coordinates already intro- 
duced this expression becomes: 


¢ 


ee ee 
(a? + yp)? 
where w is the magnitude of the volocity vector of Ps; 
with respect to P; in the rotating frame. For each posi- 


8a) it So te due (17’) 
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tion along a fixed trajectory of the Lunar Theory the 
velocity and position of the celestial body P3 will satisfy 
the above relationship with the same value of the 
constant. 

For a periodic trajectory the Jacobi integral will 
define J as a function of the period + = Qno°. The 
knowledge of J(c°) is of considerable interest in the 
study of periodic solutions and in a certain way will 
complement the equations of the orbit furnished by 
(oie 

In order to get the expression of J in terms of o we 
replace the (13) and their derivatives into (17), 
eliminate the term with fractional exponent by means 
of one of Hill’s equations and collect the constant terms. 
After some simplifications and having chosen the upper 
signs (corresponding to direct motion), we obtain 
the infinite series: 


1 =h3sd[ ene ve + ace + ve 


(18) 
9 a 9 
== 5 | Ax “ls 9 o Ay Al : 
An approximate expression up to the fourth power of 
the parameter o° can be easily obtained by using (16): 


(19) 


_ 39,533 
T76n a) 


The Jacobi integral (17’) with w = O defines the system 
of “‘zero-velocity” curves which can be used to ascertain 
the regions of space where the motion of P; can physi- 
cally occur. 

The singular points of such a system are at 


1/3 
y=0, @= +(%) 


and correspond to the value Jy = (3)*” of the Jacobi 
constant. 

For J > Jy the curves of zero velocity consist of two 
components, one of which is a closed curve T surround- 
ioe 

The interior of T is a region of stability for the 
motion of P;, meaning by that, that the third body, 
once within it, can never leave the neighborhood of the 
planet thus ensuring a permanent satellite. 

For small values of o° the Jacobi constant will take 
values larger than Jo and it is easy to ascertain that 
the periodic orbits lie within the region of stability. 


Numerical Results 


It is interesting to have numerical results from the 
series expansions already introduced. 

It is even possible to examine general characters of 
the configuration of the periodic orbits and the funda- 
mental properties of the motion of a satellite along 
them. 
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For this purpose, it is sufficient to consider second: 
order approximations for some elements like velocity 
w. radius vector R = distance (PiP3), angular dis? 
Spam 9, measured from the -AXxis of the rotating} 
system and angular velocity = @. | 

The velocity w can be obtained from the Jacobi} 
integral and R = yu *(uy) from (13): both elements} 
turn out to be infinite series of : 


cos (401°) for 1 = 0,1, 2 
T 


where r = 20° is the period of the orbit. 
Up to the second power of o we have: 


1 a Wee 3 gi ‘ 
a Ie 7 ogee 

a 6 t | ‘ 

ae cos (45 _ (20) 


(a ae yy? 


= oe (1 — : fae - *)— o COs G ‘| : 


@ can be obtained from u/v = exp (276) and (13) by 
taking the natural logarithm and considering a series 
expansion, the result is an infinite series of 


sin (40°), (T =0, 1, 2, ---), 
fy 


s 
I 


R 


I 


approximated by 


A= op ee sin (44), cl | 

T T | 

therefore, | 
. | 

oe cde peor cos (45 ) (20) 

Up ‘hi 


: See : 
The last one gives us the deviation from uniform mo-} 


The above formulae are almost self_explanatonm| 
tion, whereas the first two tell us that the velocity is) 
maximum on the x-axis and minimum on the Y-axis; | 
the radius vector reaches a maximum on the y-axis) 
and a minimum on the x-axis. These periodic trajec- | 
tories look therefore, like ovals elongated in the direc- | 
tion perpendicular to the Sun-Planet line. | 

Numerical cases have been worked out for a satellite | 
of Venus. In such case, the unit of length is the dis- 
tance Ry» = 67,270.10° statute miles; the unit of 
velocity is w*Ri. = 22 statute miles/second; the sidereal 
period of Venus has been taken as 7* = 224.7008 mean. 
solar days; and the mass ratio as w = 1/408,001. 

Evidently, the numerical value of o°, appearing in 
the formulae, expresses the period of the satellite 
around Venus as a fraction of the sidereal period of 
Venus about the Sun. 

lor the purpose of accurate numerical results, a 
program has been run on the IBM 704. The program 
essentially consists in determining the coefficients A j;x 


TABLE 1 


a3 = 0.08 J = 6.543651 

t/r % y R w @ (degrees) 
0.000 0.235928 (—2) 0.000000 0.235928 (—2) 0.300848 (—1) 0.0000 
0.025 0.232989 (—2) 0.376450 (—3) 0.236011 (—2) 0.300658 (—1) 9.1799 
0.050 0.224252 (—2) 0.743316 (—3) 0.236250 (—2) 0.300109 (—1) 18.342 
0.075 ' 0.209953 (—2) 0.109132 (—2) 0.236622 (—2) 0.299254 (—1) 27.470 
0.100 0.190474 (—2) 0.141178 (—2) 0.237090 (—2) 0.298178 (—1) 36.552 
0.125 0.166327 (—2) 0.169682 (—2) 0.237606 (—2) 0.296989 (—1) 45.581 
0.150 0.138133 (—2) 0.193961 (—2) 0.238121 (—2) 0.295802 (—1) 54.553 
0.175 0. 106602 (—2) 0.213444 (—2) 0.238584 (—2) 0.294734 (—1) 63.473 
0.200 0.725128 (—3) 0.227682 (—2) 0.238950 (—2) 0.293887 (—1) 72.348 
0.225 0.366925 (—3) 0.236354 (—2) 0.239185 (—2) 0.293345 (—1) 81.191 
0.250 0.000000 0.239266 (—2) 0.239266 (—2) 0.293159 (—1) 90.000 


numbers within parentheses denote powers of ten. 


from the recurrent procedure mentioned in the first 
part of the paper, and which is the key to any rapid 
numerical evaluation. This has been done for values of 
j from one to twelve and the corresponding values of k. 
Some of the coefficients (which we do not reproduce 
phere for lack of space) were found to be zero, so that 
while the series were summed up to the twelfth power 

of o , the maximum value of | k | appearing in the set of 

ir efiicients was six. The values of x, y, R, 6, w, and J 
were printed as functions of o° and for values of t/7 
between zero and 0.25. Good results have been ob- 
tained up to o° = 0.4. 

Table 1 is a brief sketch of numerical data for a 
satellite having a direct periodic motion around Venus. 
~¢ was chosen to be 0.08 which means that the period 

of the satellite is 17.9761 days. The unit of length is 
the mean distance of Venus to the Sun and the unit of 
velocity is the orbital velocity of Venus around the Sun, 
supposed constant; @ is measured in degrees starting 
from the x-axis, the numbers within parentheses repre- 
sent powers of ten. 


Error Analysis 

An error analysis for periodic trajectories seems to be 
essential from a practical point of view. The conditions 
which give rise to periodic orbits have been examined 
in the first part of this paper. However, it is to be ex- 
pected that trajectories which lie in a neighborhood of 
a periodic one might be of interest for the practical 
problem of tracking or forecasting the position of a 
satellite which fails to be injected into a periodic orbit 
because of engine malfunction or other causes. 

In general, these trajectories are not periodic and 
take into consideration the perturbation effects of 
higher order. 

The general theory of trajectories, lying in a neigh- 
borhood of a fixed one, can be found in [8] pages 65, 75, 
170-173, and [4] pages 335-344. Let x(t), y(t) be the 


Elements of the orbit for a satellite of Venus moving according to direct motion with period r = 17.9761 days; x, y are synodical 
oordinates with center at Venus; # is the radius vector; w is the velocity; @ the angular displacement ed from the Sun- 
Venus line; J is the Jacobi constant. Unit of length = 67 ,270.108 statute miles; unit of velocity = 22 statute miles/second. The 


coordinates of a periodic trajectory and « + Az, 
y + Ay the coordinates of a neighboring trajectory 
corresponding to the same value of the total energy 
as along the periodic orbit. The functions Ax, Ay must 
then satisfy a linear system of differential equations 
(the so-called equations of variation or Jacobi equa- 
tions) and can be visualized as isoenergetic displace- 
ments of the coordinates of the periodic trajectory. 

For every position of P3(2, y) moving along a 
periodic orbit, let us consider the intrinsic orthogonal 
reference formed by the tangent line to the orbit, 
positive in the direction of motion and the normal line, 
positive toward the interior of the trajectory. The posi- 
tion of P3(a + Av; y + Ay) moving along a neighboring 
trajectory and cor uniidite to the same value of time 
as P; can be fixed by means of its coordinates (p; q) 
with respect to the above reference. 

With the notation already used, we have 

x Yy 
peal ee 
(21) 
Ay 


i (O/ 
Deca it 


where the velocity w is obtainable from the Jacobi 
integral: 

2/3 
y) — Ju" | 


w = 2[U (a; 
(22) 


I 


bo! Go 


U(x; y) Yo + ai 


It is easy to show that the isoenergetic normal displace- 
ment g satisfies the equation 
g+ K(t)q¢ = 90, (23) 


where 


K(t) = @ ma) (1 Eft |) 
W a 


== 2 ie (Uae se Ces: 


(23a) 
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whereas the tangential displacement p is obtainable 
by numerical integration from 


p= [% (+ 


In terms of the complex conjugate variables (uw, v) we 
can write from (23a) and (2) 


noas[ied MPC) 
e/a ON 1 
Gt ; = [a are 


This expression is suitable for finding a series expression 
of K(t). Using (7’), (7), and (4) the final results can 
be put in the following form: 


: E] x = >If G cos (2 ail | (25) 


where K; (1 = 0, 1, 2, 
starting with terms in o”. 

Carrying out the calculations for the first two coeffi- 
cients it was found 


Kgs Sas 1} 
K, = 480° + Aho" + [12], 


ie dt. (24) 
w 


(23) 


K(t) 


sa 


é ; ' : 3 
--) are infinite series in a, 


(26) 


where the numbers within brackets denote powers of 
o. Neglecting the fourth power of the parameter o° we 
realize that equation (23) is of the Mathieu type. For 
sake of convenience, we let 7 = t/o° so that (23) 
becomes: 


q” +. (Ko — 2K, cos Zr) q= 0, (27) 


where the primes denote derivatives with respect to 
T. For small values of o° the coefficients Ky , K, will be 
both positive and, in general, not integers. It is to be 
understood that in order to have a complete picture of 
the situation in the neighborhood of a periodic trajec- 
tory, the general solution of (27) is needed. 

Various methods for solving Mathieu equations are 
available (see [7] page 404 and [2] chapter 5). The 
procedure followed in this paper is very convenient for 
numerical calculations when the coefficient K, takes on 
small values. 

The general solution of (27), depending upon two 
arbitrary parameters, is 


q(T; Qo; 8B) = QV Se Gi cosl (2k +c)T — Bp], (28) 


where the coefficients Q, are infinite power series in K, . 
The solution will be periodic only if ¢ is a rational num- 
ber. 

Since Ky is a small quantity, it is convenient to ex- 
press the general solution as an infinite power series in 
K, . This fact will be used to find the coefficients of (28). 

Let us remark that when K, = 0, the general solution 
of (27) is proportional to cos (1/K,7 — 8). Conse- 
quently when K; # 0 but takes on small values (this 
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is actually the situation which arises), we can write 


are functions to be determined which must necessaril 
be combinations of terms like cos [(2k + ¢)T]. 


tion and equating to zero the coefficients of the various 
powers of g, we have: c 


Bye + cE 


E,” + CE, = 2E,-1 cos (27) — F, cos (cT) 


equations with constant coefficients which will allow us 
to determine the unknown functions #;, by means of al 


q(T) = q(T; 1; 0) = cos (eT) + Ys KEE 


Ko 


I 


Ce ey Bee 
ip 
Replacing these trial series into the Mathieu equa- 


= cos [(c — 2)T] + cos [(c + 2)T] — Fi cos (cT) 


BF Ae eee ee | 
(k = 2,3, ---). (30}) 


The above written expressions are linear differential 


recurrent procedure. It is sufficient to find particular’ 
solutions of such equations; furthermore, it is easy to: 
realize that the F;’s can be determined from the} 
property that the solutions cannot contain terms like, 
T cos (cT’). If ¢ is not an integer, the above procedure 
furnishes the following result (up to the third power of 


K,): 
q(T) = cos (eT) -1K,| (c+2)T | 


c+ 1 


cos (¢c — 2)T 1 cos (ce +4)T | 
Le cancel rae ety 2) | 


| 
to=— Ie —=2) 


Lace K + eae cos (c+ 2)T 

128 (ec — 1)(e+ 1)8(c + 2) 

“ (ci 40427 Pacs(e= 2) 
(e-F Ifo) *e= 2) 
a 
3(e+ 1)(e+ 2)(e+ 3) 
COS (c—6)T 
7 ieee 
~ B(e = 1)(e = 2)(e —3) 
Roi (So erie 

Ky=e ; 
me +3@-1) * Bee —D 
The second relationship, written above, should be 
used to determine the value of the parameter c; how- 
ever, it is not very convenient for numerical colonial 
tions since a high order algebraic equation in ¢ should 
be solved in order to get a good approximation. We 
shall follow a different procedure to obtain the value of 
c explicitly in terms of Ko and K,. From the general 
theory of Mathieu equations, it is known that if 


+ 


(31) 


+... 


p(T; Ko, Kx) is the particular solution of (27), such 
hat ; 


mG;-Ko, Ky) = 1 ~and Cite 1 CEG 0, 


COS (ec) — g(a: Ko ) ky), (32) 


.e [1] page 208 and [3] page 103. 
Expanding g* into a power series of K, 


q*(P; Ko, Ki) = Di, Ki'g(T; Ko), (33) 


Delta GT ie, oi) 
r be K — 2 2 
q ( 0) r! | OK" vee 


go(T; Ko) = cos ~/KyT, (34) 


whereas the other q,’s satisfy the recurrent. system of 


, (r = 1, 2, 35 ns (35) 
with initial conditions 


‘ 
. qr(0; Ko) = qr'(0; Ko) = 0. (35") 


Hence, up to the fourth power of kK, it was obtained 
from (32), (34), and (35): 


cos Pet cos(r4/K,) + Ky (rV/Ko) 


4./Ko(Ko — 1) 


. KS [ee — 32K, + 8)msin (1+/Ky) 
: 64(Ko — 1)?(Ky — 4) Ko VK, 


- rT COS — (36) 
32K (Ko — 1)? | 


This formula gives explicitly the value of cos (ze). 
The values of c thus obtained differ by an arbitrary 
additive constant +2k (where k is an integer); the 
value to be chosen must be the closest to ~/Ky. In 
general, c is not an integer. If this is not the case, the 
solution of the Mathieu equation (27) can be obtained 
from standard procedures available in many textbooks, 
see e.g., [2] pages 15-17 and [7] page 404. 

The results expressed by formulae (31) and (36) 
complete the required elements for the general solution 
(28) of (27), the two parameters Qo and 6 being deter- 
mined from the initial conditions (i.e., position and 
velocity ). The knowledge of q(t/ o ) and of the elements 
of a periodic orbit (x, y, w and their derivatives) allow 
us to obtain the tangential displacement p from (24) 
by means of a numerical quadrature. Equations (21) 
finally yield Ax and Ay in terms of the time. 

Equation (23) is the well-known Hill equation and 
turns out to be of fundamental importance in the 
problem of the stability of periodic orbits, a problem 
which is strictly connected with the subject matter of 
his section. Among the numerous literature on this 


topic one should mention a recent paper yan ame 
Message “Some periodic orbits in the restricted prob- 
lem of three bodies and their stabilities” published in 
the Astronomical Journal 64, 226, 1959. 


Conclusions 


The trajectory of a stable satellite around a planet, 
when the Sun’s perturbation is taken into considera- 
tion, should be a periodic one. 

The periodic trajectories obtained in this paper from 
the point of view of the Lunar Theory, are symmetrical 
with respect to the axes of a synodic reference with the 
origin at the planet; they are oval-shaped with the 
planet at the center of symmetry and slightly elongated 
in the direction perpendicular to the Sun-Planet line. 
The coordinates of such trajectories have been ex- 
pressed as cosine and sine series of time. A recurrent 
procedure allows us to obtain the coefficients of the 
expansions, once the period of revolution has been 
chosen, with any required degree of accuracy; the 
procedure is furthermore fit for programming on a 
high-speed electronic computer. The Jacobi constant 
of the motion can likewise be computed in terms of the 
period, so that from the Jacobi integral it is possible to 
obtain the magnitude of the velocity for positions along 
an orbit of such period. The numerical value of the 
Jacobi constant is therefore indispensable for evaluating 
the requirements of a retrorocket so that a ballistic 
missile can be injected into a periodic orbit. 

The numerical work undertaken for a satellite of 
Venus yields good results for trajectories whose period 
is up to 0.4 of the sidereal period of Venus. To account 
for any practical situation, deviations from a periodic 
orbit have been studied as displacements of the posi- 
tions of the periodic orbit itself. The numerical methods 
for such calculations expressed in the paper seem quite 
straightforward for applications. 
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On Stromgren’s Method of Special 


Perturbations 


Peter Musen* 


Abstract 


In this paper Strémgren’s method for the numerical com- 
putation of perturbations is modified to include the effects of 
higher order. The accurate form for the rotation matrix in 
terms of, so called, Gibbs vector is given and the differential 
equation for the perturbations of this vector is established. A 
similar modification is applied also to Hansen’s method. 


Introduction 


In this paper Strémgren’s method (1929) for com- 
putation of special perturbations is modified to include 
the effects of higher order. The standard form of Strém- 
gren’s method possesses a mathematical elegance and 
simplicity and it gives the vectorial elements of the 
orbit without trigonometrical transformations. The 
disadvantage of this method is that it takes into account 
perturbations of the first order only. The integrated 
value of the angular velocity of rotation of the osculat- 
ing ellipse is used by Strémgren to obtain the Gibbs 
rotation vector, and because this vector is not ob- 
tained directly, only the first approximation to the 
matrix of rotation can be easily deduced. 

The author decided to reconsider the possibility of 
returning to Strémgren’s method and to deduce the 
accurate equations for the Gibbs vector. Using this 
vector a closed form for the rotation matrix is also 
obtained. 

Some methods developed by S. Herrick (1948) and 
by the author (1954) avoid the forming of the matrix 
of rotation by using three unit vectors rigidly con- 
nected with the osculating ellipse, but, in these methods 
the eccentricity appears as a ‘small divisor’ in more 
equations than in Strémgren’s method. This fact 
limits the applicability of all three methods to case of 
moderate or large eccentricities. The new form of 
Herrick’s method (1953, 1961) avoids the division 
by e completely and can be used also for the cases of 
small eccentricities. 

The division by e does not appear also in Hansen’s 
method. The form of Hansen’s method which is con- 
venient for the computation of special perturbations 
for all eccentricities is given in the last part of this 
article. 


* Theoretical Division, Goddard Space Flight Center, 


National Aeronautics and Space Administration, Washing- 
ton, D. C. 
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The Basic Equations 4 
Let R be the unit vector normal to the osculat 
orbit plane, P be the unit vector directed from 1 
central body to the osculating perihelion, and Q 
R <P. The values of these vectors for the ini 
moment will be designated by Po, Qo, Ro. As in = 
author’s previous article (1954) the equation of 
motion is taken in the form 
dr r 
dé 7 E, 


. 
and the projections of the disturbing force F on | 
Q, R, rand R X rare designated by L, N, W, S and 
respectively. Let 3, ¢ and wy be the inclination of t 
orbit, the longitude of the ascending node and t 
argument of perigee of the osculating orbit, where t 
system Po , Qo, Ro is taken as a reference system. ; 


cost (g¢—y), 1 | 
ieee ig 8 2 
"cost (ety /2 j 

sins (¢ =) 
= ig 08 
7 cos (oy) 92 
gos=tgz(e+y) (- 
gs = mPo + geQo + gsRo. € 


The matrix of rotation T, transforming the syste 
(Po, Qo, Ro) into the system (P, Q, R), can be wri 
ten in the forms suggested by Gibbs (1901) 


or 


1-¢” 2g X 1+ 2g¢ 
fi r+ : 
1+ ¢ l+¢ ; 


where J is the unit matrix, 
IT = PoPo + QiQo + RoRo. 


The equation (6) or (6’) replaces the development 
into an infinite series (Herget, 1948). Let e@ and 
be two Gibbs vectors. If the rotation defined by y fi 
lows the rotation defined by g, then Gibbs vector 
the combined rotation is given by the formula (GibI 
1901). 


(6 


‘ps See eS 


6 
Li iguy ( 


Phe osculating unit vectors P, Q, R are given by the 


Za aan eo Us eX (2 XU) 
tks 1—¢ 2e X U) + 2g¢-U, 7) 
Up s 0 ss 0 


1+ ¢ 1+¢ : 


where P, Q or R must be substituted for U and Py, 
Qo, Ro for Up. After P, Q and R are obtained the os- 
culating position vector is computed using the stand- 
ard formula: 


r = Pa(cos EH —e) + QavV/fl—e@sinE. (8) 


Let Q be the instantaneous angular velocity of rota- 
tion of the osculating ellipse and w; , w: , w; the projec- 
tions of Q in the directions of P, Q, R, we have (Strém- 
Seren, 1929): 


w, = hWa(cos FE — e) (9) 
w =hWarv/i—ésinE (10) 
, w= V4 (sing - V1 —éL) ‘Ga8) 
4 
J — J — , 
: V/a(l — &) 
_ Let m be the unit vector along the line of nodes in the 
: system (Po, Qo, Ro). It is kinematically evident that 
er oO + oR 
_ dd 1 dle +¥) 
1 d(g — p) 
a 5 (R, — R) Soret 
Put 
@ = @1Po + Qo + wsRo- (13) 
We have 
Q2=T-e=ot+ jlgXo+gX (gXo)]. (14) 


=a 


Let g and g + dg be two Gibbs vectors defining the 

rotations of the osculating ellipse from the moment 

fy to the moments ¢ and ¢ + dt respectively. The 

infinitesimal matrix of rotation which brings the os- 

culating ellipse from the position at the moment ¢ 
to the position at the moment ¢ + dt is 


IT+od XI, 


and the corresponding Gibbs vector is 4 dt. Substitut- 
ing 


G=e+dg, g=8- Y¥ = 70d, 
into (6”), we deduce 
2 ho Age i ees : 
. Fae Og See ty BETO (15) 


OF 
dg _ 
dt 


We find from (6’) 
aig Te Pee er eee 
From (17), (16) and (14) we deduce 


s(I+eXI+e8)-0. (16) 


dg _1 2 

aq et g(a + o), (18) 
and also 
S(t @e tse Xxotse x (exe). (19) 
dt g +58 0) ee g Xo 


Any one of the equations (15), (18) or (19) gives 
Gibbs vector g accurately and replaces Strémgren’s 
equations for dé,/dt, d6,/dt, edr/dt. If the second and 
the third terms on the right hand side of (15) are 
neglected, the result coincides with the basic Strém- 
gren’s equation. Multiplying (16) by the matrix 


Ast I 
eee Seay 


and considering the following formulae, 
I-A=A 
ALXEISR WAG: 
(g X I)-(g XI) = gg — gl 


we deduce 


Eye 2 dg 
Cae baege ed Jones 
or 
__2 (dg, dg 
al: i a x), (20) 


This vectorial formula is identical with the classic 
scalar formulae (Whittaker, 1904): 


E dy dé dx 

a =2(+x 84 pit _ _) 
Pay Fibra. dn dg dx 
a= 2 a eae ee 0 tx) 


dt dx 
a, = 2 = we 
2 ( pe dt eit TX dt s 


where & n, ¢, x are Euler’s parameters 


~~ 


i bre Gin 
é = sin 5 3 cos 5 beer ee et 
Vit+¢ 

n= sin 5 : 5 OSD = 5 ( —y)= ae LS 
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1 
(G = Cals sin 5 
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ie: 
(go+y) = 
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We have for the perturbations of the mean anomaly 
(Brown, 1896) 


dAM 2rS ——— 
COE agen Se Aenea (21) 
dt Ve ieee ae 


Multiplying (12) by Ro + R, we deduce 


d(g == y) - (w,P + wQ)-Ro(1 at g) (22) 
i ee 2(1 + gs) 

or, taking (7) into consideration, 

dge+y) _ pees w1(Gigs — 92) + w2(gegs + 9) (93) 


near 1+ gs 
Adding (21) and (23) together we obtain an equation 
for the mean longitude ZL; in the system (Po, Qo, Ro), 


aly ee 2r8 
dt ve ete Te ver = =; (ee) 
4 (Gags = go) + wlgags + 91) 
ios gs 


This last equation replaces Strémgren’s equation for 
dL,/dt. The equations 


(24) 


@ = Vall —@) (Leos E +N) (25) 
vi EN Aho ARG EE ot (26) 


dt avV/1— é 
remain unchanged. Kepler’s equation takes the form 
E —esinE = M,+ (AL, — 2 arctg gz) 


2g 
+ no (t= 0) + ff Fae, ae 


The values (9) — (11) of 1, we, w3 must be used in 
the equations (15) — (24). The equations (15), or 
(18), or (19) and (24), (25), (26) and (27) constitute 
the basic system to be programmed for the use of elec- 
tronic machine. In this modified Strémgren’s method 
the danger of losing the accuracy because of the small 
divisor is not greater than in the classic Lagrange 
method and the matrix of rotation is formed directly, 
using three independent parameters. The developed 
equations can be easily programmed and the author’s 
opinion is, that it is desirable that Strémgren’s method 
come into wider use again. 

A similar system of formulae can also be developed 
for the use of Hansen’s method (1857) to compute the 
special perturbations. Let P, Q, R now be the system 
rigidly connected with the osculating orbit plane (not 
with the osculating ellipse as in Strémgren’s method). 
The basic Hansen’s ellipse has the fixed value ay of the 
semi-major axis and of the eccentricity e) and it does 
not move relative to the system (P, Q, R). Let 1 + » 
and AM be Hansen’s perturbations in the radius-vector 
and the mean anomaly. The position vector r of the 
disturbed body is given by the equations 


= Ce p (28) 
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= Pacos EF — é) + Qa l—ezsni (29 

E —esin E = Mo + no(t — t) + AM. (30) 

The instantaneous angular velocity of rotation of the 
system (P, Q, R) is | 
= (1 + v)hWo. ’ (31) 

Put 
o = (14+ ~)hWo, - (32) 
= Pya(cos H — eo) + Ove —ezsin H. (33) 


ie ieee 
AW 
dt 4 Co + 


and evidently, 


+ o)(1 ac v) (34) 


= 00 a lg X o +8 X (g X oo). (35) 


meg 


The remaining system of formulae can be retained 1 
the form given in the author’s previous article (1958) 


w= X + Y (cos E — @&) 


(36) 
+ Z/1 — e? sin E | 
ho 1 
h ; V/ao(1 a eo”) 
dA = Al 4 pT (38 
dt 
oie fe C1 + y)- T cos Bi x] ; 
dt he? 
dZ ae ey ore | 
=e Oh: . vt E-L (40 
= hy A E Noe we sin 
Z (8+ A)A 41 
“ En (41) 
1 10 AN Y 
ts eee} =——s 42, 
1l+p r 2(1 + A) : 


and (Hill, 1881) 


d\M_ wty 


— hy 43° 
dt a 1— pv (43, 


The system of formulae (34) — (48) does not contair 
the eccentricity as a divisor and, consequently, it car 
be used for all orbits. This system was programmed a‘ 
the Minor Planets Center in Cincinnati for the com 
putation of special perturbations of celestial bodies 


Acknowledgement 


The author takes this opportunity to express hi 
gratitude to his colleague Mrs. Ann Bailie for carefu 
checking of the manuscript. 


References 


1. Brown, E. W., Introductory Treatise on the Lunar Theory 
1896, pp. 58-63. 


2. Gipss, J. W., Vector Analysis, Yale Univ. Press 1901, 
pp. 343-345. 

3. Hansen, P. A., “‘Auseinandersetzung” I, Abh. Sach. Ges. 
Oicnad) Wiss. 5, 1857. 

4. Hercer, P., The Computation of Orbits, 1948, pp. 84-85. 
>. Herrick, S., Pub. A. S. P. 60, 1948, pp. 321. 

§. Herrick, S., ‘Icarus and Variation of Parameters,” A. J 
58, 1953, pp. 156-164. 


hd 


Abstract 


This paper traces the career of Dr. Robert Hutchings 
roddard from his earliest rocket experiments as an under- 
graduate at Worcester Polytechnic Institute to his work at 
Jamp Devens, Massachusetts. It also attempts to demon- 
trate that during this period Goddard was recognized by his 
eers, that he was supplied with funds for research and de- 
‘lopment, and that he was a very thorough and careful 
eientist who eschewed the sensational in his work and find- 


ngs. 


ntroduction 


This paper is almost in the nature of an attack, an 
tack on certain near legends and myths regarding the 
areer and achievements of Dr. Robert Hutchings 
xoddard, late Professor of Physics, Clark University, 
Worcester, Mass. Over the years since his death at 
3altimore, Maryland, there have grown stories to the 
ffect that he was ignored by his peers, that he was 
aughed at by the knowledgeable, that he was without 
unds, that he labored alone and unknown, that he was 
/ prophet without honor in his own country. 

Like so many legends and myths there is a grain of 
ruth in these stories, but for the most part they simply 
re not true. If this paper accomplishes nothing else, 
f will show that Dr. Goddard was well known in his 
ime; that he was amply and very often supplied with 
esearch grants, funds, and facilities from the Smith- 
onian Institution, the National Geographic Society, 
he Carnegie Institution, Clark University, Worcester 
*olytechnic Institute, and especially the Daniel and 
‘lorence Guggenheim Memorial Fund; that he did not 
ibor alone; and that the public and the scientific com- 
iunity were almost always well aware of what he was 
bout (the New York Times never failed to publicize 
im, often on Page One; as a matter of fact, his widow 
as in her possession some nine volumes of newspaper 
lippings on her husband). 

It is hoped, too, that this paper accomplishes some- 
ng in the way of destroying a pejorative tale about 
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Goddard and His Karly Rockets: 
1882—1930 


KE. R. Hagemann 


Goddard: that he was a “moony” man; that is, he 
emphasized the extravagantly romantic and somewhat 
unscientific aspects of rocketry, that he was a Sunday 
supplement writer in the vein of a present-day scientist 
who need not be named. Goddard, except for a minor 
lapse or two, was precisely the opposite. He was a very 
careful, very thorough scientist who although suggest- 
ing an ultimate result—a landing on the Moon—was 
always concerned with goals nearer at hand. He was, 
in his own way, a systems engineer. He was an accurate 
and careful mathematician. Unfortunately, the way- 
ward press of the United States inevitably extracted 
the sensational from his work and stressed it. 

This horrified Goddard. It stiffened him. It made him 
suspicious. He refrained from publishing (only three 
articles on rocketry were published in his lifetime). But 
he never slackened in his tasks and he worked steadily 
until his death, as he had worked all his life, on the data 
at hand. 

Although he appears to have been a pragmatic, prac- 
tical man, he was not without a romantic aura, par- 
ticularly in his early life; and this places him in the tra- 
dition of the Yankee inventor, Edison, Goodyear, Fire- 
stone, Ford, et al. He was born in Worcester, Massa- 
chusetts, on October 5, 1882. Because he was a sickly 
and thin child, he was prohibited from the usual aimless 
violence of a boy’s life. Instead, he read voraciously and 
nursed his imagination. In his middle age, Goddard re- 
called at the age of six having seen electricity produced 
by a Leclanché battery. He obtained the zine from such 
a battery and scuffed along a gravel walk. He mounted 
a low fence and jumped. Then he repeated the experi- 
ment ‘‘scuffing over a longer distance, and endeavored 
to convince myself I had jumped higher.” 

In his sixteenth year, he says, “‘an event happened 
which was destined to provide me with all the scientific 
speculative material I could desire.” The Boston Post 
in January 1898 serialized Fighters from Mars; or, The 
War of the Worlds, In and Near Boston, and followed 
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with Garrett P. Serviss’ Hdison’s Conquest of Mars. He 
also read Jules Verne and was very fond of From the 
Earth to the Moon (1865) in which rockets were featured. 

His teen-age mind roamed and probed without rest. 
In the Spring of 1898 he thought it would be perfectly 
delightful to fly a permanent balloon at the end of a 
string. He obtained 1/100-inch aluminum sheet, shaped 
it like a pillow, and filled it with hydrogen. It was natu- 
rally too heavy. He accepted his “failior,” as he spelled 
it, and turned to other things: “how birds fly,” for ex- 
ample. On October 19, 1899, he sat in a tree and “im- 
agined how wonderful it would be to make some device 
which had even the possibility of ascending to Mars.” 
His crude wooden models gave “negative results,” and 
Goddard admitted ‘that there might be something 
after all to Newton’s Laws, which I had read in Cassell’s 
Popular Educator, given to me by my father.” 

He was late in entering high school but not inactive 
there. During his first year at South High in Worcester, 
in 1901, he wrote an article on space navigation and 
submitted it unsuccessfully to Popular Sccence Monthly. 
Goddard writes: ‘The article mentioned how the idea 
of firing several cannons, arranged like a nest of beakers, 
had been proposed, and quoted figures on the frequency 
of meteors in space... .”’ It will be noted how he always 
considered various possibilities and problems. (I have 
mentioned his thoroughness and his carefulness.) In 
other words, he measured. Propulsion by firearms of 
some kind was a favorite idea of his, and he once con- 
sidered obtaining propulsion by a machine gun device 
that fired bullets downward—he cannot forget Newton’s 
third law. 

Upon graduation from high school in 1904, Goddard 
heroically burned all his notebooks because “I knew 
they were erroneous.”’ But the dream would not be for- 
gotten and inside of two months he was again scribbling 
and calculating. He bought a batch of green clothbound 
notebooks and commenced a scientific and systematic 
recording of his ideas and concepts. (They are fasci- 
nating and we will return to them in a moment.) 

Goddard matriculated at Worcester Polytechnic 
Institute and was graduated from there with the degree 
of Bachelor Science in 1908, holding highest honors in 
physics and mathematics. His academic record did not 
hold him in thrall: he was president of his class, editor- 
in-chief of the college yearbook, and a fraternity man. 
Nor did the textbooks hold his imagination and curi- 
osity in thrall. In his freshman year a professor assigned 
“Travelling in 1950” as a theme topic. Goddard re- 
sponded with a bold and brilliant idea: a railway line, 
with air-tight cars in a vacuum which were supported 
electromagnetically without metal-to-metal contact. He 
calculated that there would be an acceleration of 32 ft/ 
sec? during the first half of the journey and a decelera- 
tion of equal amount during the last half. This meant a 
running time of ten minutes from Boston to New York 


16 vols. Ist pub. 1854. 
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cerned with the demonstration of a principle: the | 

tinued acceleration of a body by forces which chan 

from attraction to repulsion as the body passed by } 
source of the force. 

It would seem that propulsion was almost an obs 
sion with him. In his junior year at WPI he submitts 
seriatim, a paper to Scientific American, Popular © 
ence Monthly, and Popular Astronomy that propoy 
heat from radioactive materials be-used to expel 8 
stances at high velocity from an orifice. Jet propuls 
would be furnished sufficient to navigate in interplag 
tary space. The periodicals refused the manuscripth 

Toward the end of his senior year, beneath a build 
at Worcester Tech, occurred the most significant eve 
of all: a smoke-filled basement, the result of static t i 
of small rockets. | 

I now must return to his notebooks. They are exte 
and furnish us a classic example of usage by a Et! 
grounded man. They also illustrate how a notebe 
might fecundate a mind. 

During the Christmas recess in 1909, Goddard st 
marized twenty-six methods in his notebooks “4 
ing means in space, means taken with the apparat| 
and means sent from the Earth.” There isn’t the ro 
for all of them here; however, extracts suggest thee 
ness thereof: 

(1) September 6, 1906: reaction by streams of ions# 
furnish rocket propulsion. (This was subj 
quently developed in a patent taken out May 
1915.) 

(2) July 1907: use of solar energy in connection y 
electrostatic repulsion. (Goddard was to do f} 


City. Typically enough, however, Goddard was 


Goi sae ) 

(3) April 4, 1908: large mass of explosive to raise ( 
pound final mass to great heights; with slow pI 
pulsion in the atmosphere. (This is, of count 
germane to his work in 1919.) . 

(4) June 9, 1908: liquid hydrogen, oxygen, nite 
feeonte (N2O,), and ethane (CsHe). 

(5) June 19, 1908: camera sent around distant plan 
and returned to Earth. 

(6) June 24, 1908: circling a planet to decrease | 
wats landing. 

(7) October 15, 1908: steering automatically by pho! 
sensitive cells. 

(8) January 24, 1909: multiple rockets. (He was to. 
further work on this while honorary fellow | 
Clark.) 

(9) April 6, 1909: cooling a nozzle by liquid hydrog 
and oxygen. 

Within a year, Goddard had added such ideas as p 

duction of hydrogen and oxygen on the Moon, the ge 

eral theory of the hydrogen and oxygen rocket, and 

means for neutralizing the effect of the decrease 

gravity on an operator. 

The six to seven years subsequent to graduation fre 


| 
| 


PI, that is to 1915, were intense and trying for 
ddard. He taught and he continued graduate work. 
obtained an A.M. from Clark University in 1910 
d his Ph.D. in physics from the same school in 1911. 
‘doctoral dissertation was the explanation of crystal 
ifiers; or, to use its formal academically heavy title: 
n the Conduction of Electricity at Contacts of Dis- 
nilar Solids.” (It was published in the Physical Review 
June 1912.) Goddard admitted candidly that his 
toral subject was chosen not because he was inter- 
ed in it. Far from it, apparently. But because his 
vious studies ‘‘on the conductivity of powders” at 
T would help him. He squandered nothing! 
3y now Goddard was combining teaching with doc- 
research and postdoctoral research. He was in- 
ctor of physics at WPI from 1909-1911. He then 
ent a year at Clark, 1911-1912, as an honorary fellow 
physics. From 1912-13, he was a research instructor 
physics in Princeton University. In 1914, he returned 
y Clark, where he was to remain academically, as in- 
ructor and fellow in physics. In 1915 he was promoted 
stant professor and 1919 promoted professor of 
sics. In eight years, not all of which he spent in resi- 
¢e, he had advanced from Ph.D. to full professor. 
These very early years as teacher and research fellow 
7 important to our understanding of Goddard’s total 
reer. When at Clark in 1911-1912, he worked on vari- 
is rocketry methods, and the theory of the multiple- 
large rocket was roughly outlined as was the use of a 
rhtweight solar energy engine. (These were mentioned 
his notebooks.) At Princeton, he was engaged day 
id night. He writes: ‘During the day I worked on an 
teresting physical problem not then explained by 
setrical theory—the positive result of force on a ma- 
rial di-electric carrying a displacement current.” This 
ight be explained a little more familiarly by saying he 
oduced the first laboratory demonstration of me- 
anical force from “a displacement current” in a mag- 
tic field. This is a fundamental concept in Maxwell’s 
eory of electromagnetic waves. In the evenings, 
sddard worked on rocket propulsion theory. He as- 
med that an efficiency of 50 per cent could be secured 
th nitrocellulose smokeless powder and hydrogen and 
ygen. This theoretical work was the basis for the 1919 
tblication. 
Byen when the schedule became too much and he was 
Miously ill with tuberculosis in 1913-14, Goddard did 
i really rest and in May 1913 he wrote the material 
*two US Patents which he himself later said covered 
2 essentials of rocket propulsion. They were taken out 
thin a week of each other in July 1914: No. 1,102,653 
Rocket Apparatus (2-step Powder Rocket) and No. 
103,503—Rocket Apparatus (Cartridge Reloading). It 
better that we allow Dr Goddard to speak of them: 
. They give us as nearly as possible an answer to 
he question as to what the “Goddard Rocket” 1 


..l prefer...to consider the three niceiples 


covered in the claims of these two patents, namely 

(1) the use of a combustion chamber and nozzle; (2) 

the feeding of successive portions of propellant, liquid 

or solid, into the combustion chamber, giving either 

a steady or continuous propulsive force; and (3) the 

use of multiple rockets, each discarded in succession as 

the propellant it contains is exhausted. 

In the Fall of 1914 Goddard returned to Clark as a 
part-time teacher and fellow and immediately began 
further research that would, within five years, lead to 
his now famous 1919 Smithsonian publication. Al- 
though his entire career to 1919 was a preparation, the 
period 1914-1919 was particularly so. That Fall he 
worked out completely the theory and calculations for 
smokeless powder, hydrogen, and oxygen and began ex- 
periments on the efficiency of ordinary rockets. He 
learned that the initial mass needed to send one pound 
to infinity for hydrogen and oxygen, at 50 per cent effi- 
ciency, was 43.5 pounds, very close to the 45 pounds 
estimated on January 31, 1909. 

The pace of his work at Clark was stepped up. He 
experimented with a vacuum tube “in which ions 
moved in closed paths due to a magnetic field.”” On 
May 4, 1915, he obtained a patent for ‘““A Method of 
and Means for Producing Electrically Charged Par- 
ticles.’”’ In September he patented an apparatus for 
producing gases, hydrogen and oxygen, for manufac- 
turing liquid oxygen and hydrogen. 

In 1915, Goddard began building solid propellant 
rockets and over the next couple of years was to spend 
some $800 of his own money to support these experi- 
ments. Chiefly, he was concerned with measuring the 
efficiency of common rockets, the efficiency of steel 
rockets provided with nozzles, and efficiency of steel 
rockets in vacuo. With great care, proceeding very 
cautiously, this work was written up. By the middle of 
1916, he had reached the limit of what he could ac- 
complish on his own resources. He began to look around 
for funds. 

He gathered together some photographs with his 
notes, had them bound in a book with a green cover 
and a neat gold border, and called the manuscript, 
“A Method of Reaching Extreme Altitudes.” Its aim 
was simple: to secure sufficient funds to permit the 
development and construction of a rocket with a large 
proportion of propellant to total weight. 

It might be educational to observe Goddard’s cau- 
tious exclusions and inclusions in the manuscript. For 
example, although measurements for projection to in- 
finity were included, to avoid giving the impression of 
too much speculation, he thought it best to suggest 
sending a mass to the Moon rather than suggest inter- 
planetary space navigation. He omitted the discussion 
of solar energy to decrease transit between planets. He 
omitted the discussion of the Moon as a half-way sta- 
tion. But he left just one tantalizing hint: 

There are, however, developments of the general 

method under discussion which involve a number 
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of important features not herein mentioned, which 

could lead to results of much scientific interest. 

A letter apropos the manuscript was dispatched to 
the Aero Club of America. No reply. Other letters 
followed. No encouragement. Then, almost at the end 
of his list, Goddard sent a letter to the Smithsonian 
Institution. He received an encouraging reply and a 
request for the manuscript to be examined by a com- 
mittee. It was December 1916. 

In about three weeks, Goddard received a letter from 
Dr. Charles D. Walcott, then Secretary of the Smith- 
sonian, commending him on the manuscript and asking 
him how much money was needed. With typical aca- 
demic caution and wariness, Goddard debated between 
asking for $2500 and $10,000. He settled on $5,000. 
In the next letter from the Institution there was an 
advance, a check for $1,000, more money that he had 
ever seen at one time. 

Then began early in 1917 the series of experiments 
which would, we must know, launch seriously modern 
rocketry and which would bring fame to Goddard. The 
work was carried on chiefly in the Magnetic Laboratory 
at Worcester Polytechnic Institute but, also, in part at 
Clark. Such work had no more than commenced, when 
World War I was declared. Goddard desired to under- 
take the work as a defense measure and in the winter 
of 1917 he obtained approval in Washington and was 
soon named Director of Research, U. 8. Signal Corps. 
Until June 1918, research continued at WPI; from 
June until almost November it was carried on at the 
Mount Wilson Observatory shops in Pasadena, Cali- 
fornia. 

On November 10, 1918, at the Aberdeen Proving 
Grounds, before representatives of the Signal Corps, 
U.S. Flying Corps, Army Ordnance, et al., Goddard 
demonstrated the results of his research. He launched 
an anti-tank and anti-personnel projectile rocket, 
about three inches in diameter, from a launching tube 
with two short rear legs. It was the forerunner of the 
bozooka of World War II. He also launched a trajectory 
rocket, a small multiple-charge rocket, with five charges 
in a combustion chamber much like a repeating rifle. 
It carried a simulated warhead. It travelled “straight,” 
Goddard says, for 60 to 80 feet. 

The Armistice next day put an end to Goddard’s 
work for the Army. 

Upon his return to Clark in 1919, Goddard was 
urged by a colleague to published his data. He turned 
to the Smithsonian, asked them. He correctly thought 
a single publication more effective than several papers— 
an instructive idea. The Smithsonian agreed to publi- 
cation zf cost were taken from his grant for rocket 
development. Goddard assented. Not quite as cautious 
as he was with the 1916 manuscript, he included supple- 
mentary discussions on secondary rockets, use of 
hydrogen and oxygen, meteoric collisions, ete. 

Full bibliographic information is as follows: 

Goddard, Robert H. “A Method of Reaching Ex- 
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treme Altitudes.” Smithsonian M iscellaneous C 
lections, 71, No. 2 (May 1919), 69 pp., with) 
plates. Publication Number 2540. | 

By now it is a collector’s item. Exactly 1,750 cops 
were published. Although dated May 26, 1919, 1t Ay 
issued simultaneously with a news release by the Ing 
tution on January 11, 1920. It is basically the saq 
report as that of 1916. Goddard’s experiments sim} 
corroborated his previous conclusions; only the fac 
data based on his post-1916 tests needed revision. If 
work had cost a total of $11,000. The paper was divig 
into three main parts: Theory, Experiments, and Cala 
lations Based on Theory and Experiment. There w i 
seven appendices containing supplementary discussio | 
additional notes, and 29 photographs. q 
Early in the paper, Goddard stresses the importa 
of the subject and suggests a means to explore tf 
atmosphere above 20 miles. ‘In fact,” he says, ™ 
most interesting, and in some ways the most importai 
part of the atmosphere lies in this unexplored — 


to be investigated were density, chemical constitutic 
temperature, the aurora, and alpha, beta, and gama 
rays. He then makes this curious statement: “ 
instruments for obtaining data at these high altituc 
are herein discussed, but it will be at once evident th 
their construction is a problem of small difficulty CoH 
pared with the attainment of the desired altitudes.’| 

We should accent that Goddard’s prime interest W 
this region 20 to 50 miles above Earth; we should 
cent that its exploration and its scientific analysis w r 
uppermost in his mind and not planetary or int} 
planetary travel. | 

Regarding rocket action, Goddard continues, “4h 
problem was to determine the minimum initial mass: 
an ideal rocket necessary, in order that, on continug 
loss of mass, a final mass of one pound would remd 
at any desired altitude.’’ Provided that a velocity} 
7000 ft/sec could be attained and that most of t 
rocket was propellant, to raise one final pound to | 


miles required a starting weight of 3.66 pounds; | 
send one pound to infinity, 602 pounds. | 

To hit the Moon (and with this obiter dictum Gc 
dard titillated the newspapers), a total initial mass) 
8 or 10 tons ‘‘would, without doubt, raise sufficie 
flash powder for clear visibility.” It was his plan 
send Victor flash powder to signal impact with t 
lunar surface. 

After initial experiments with black powder, | 
used dense smokeless powder. Charges were fired 
steel chambers with throat diameters of about or 
quarter and one-eighth inch and 8-degree diverge 
nozzles with various expansion ratios. An efficiency 
over 64 per cent and an average velocity slightly unc 
8,000 ft/sec were attained. Then to prove these 4 
locities real and knock down a prevalent unscienti 
notion, Goddard fired the same steel chambers in vac 
and observed the recoil. There was a 12 to 20 perce 


mcrease in velocity. He observed: “If... successive 
harges were fired in the same chamber, much as in a 
id fire gun, most of the mass of the rocket could con- 
sist of propellant... .” This was what Goddard termed 
“step-rocket”’ and discussed it to some extent. 

He drew an early conclusion: “It is believed that 
ot only has a new and valuable method of reaching 
igh altitudes been shown to be operative in theory, 
4 the experiments herein described settle all the points 


Although one critic in the 1950’s said that “A 
Method” “did not indicate a talent moving with ease 
in the field of mathematically controlled speculation,” 
Mr. Alfred Africano (STL), who subjected the paper 
to a thorough analysis in 1946 and applied Goddard’s 
¢alculations to the V-2, found an impressive ‘“‘agree- 
ment of the actual with the calculated minimum weight 
ratio within 2 per cent for the first interval [5,000 ft], 
1 per cent for the second [15,000 ft], and 3 per cent for 
the third [25,000] and fourth [45,000].”’ 
~ (Intervals were Goddard’s way of dividing altitude.) 
_ Too bad, really, that Goddard was not removed to 
one of his intervals when the story of Publication 2540 

roke on Page One of the New York Times, Monday, 

fanuary 12, 1920. 


BELIEVES ROCKET 
CAN REACH MOON 


This was the headline. The story related how God- 
dard had invented a ‘‘new type of multiple-charge, 
high efficiency rocket of entirely new design”’ to explore 
“the unknown regions of the upper air’ and even the 
Moon itself. Then followed a lengthy summary of the 
paper in which the scientific elements were stressed, 
something the headline did not do. The Times felt 
that sending recording apparatus to moderate and 
extreme altitudes would greatly benefit meteorology 
and weather forecasting. Indeed, the story said boldly 
that a Moon shot ‘“‘would be of little obvious scientific 
value,” although of “great general interest.” 

The Olympian Times had a good thing. In ‘“‘Topies of 
the Times,” January 13, the editorial writer, under the 
sub-head, ‘‘A Severe Strain on Credulity,” was nasty: 

. After the rocket quits our air and really starts on 
its longer journey, its flight would be neither ac- 
eelerated nor maintained by the explosion of the 
charges it then might have left. To claim that it 
would be is to deny a fundamental law of dynamics, 
and only Dr. Einstein and his chosen dozen, so few 
and fit, are licensed to do that. 

The tirade delivers itself of fustian: 

That Professor Goddard, with his ‘‘chair” in 

Clark College...does not know the relation of 


B 
a 


, 


action to reaction, and of the need to have something 

better than a vacuum against which to react—to say 

that would be absurd. Of course he only seems to 
lack the knowledge ladled out daily in high schools. 

A newspaper should know when it has a story. On 
January 16, the Times printed a letter from crusty 
Admiral William 8. Sims, USN, who insisted that a 
rocket would ascend in a vacuum. He recalled a physics 
experiment at the Academy in 1882. 

Goddard himself was aghast and wry. In the public’s 
mind, it was a ““Moon-rocket.’’ And Goddard felt that 
“by trying to minimize the sensational side, I had 
really made more of a stir than I would if I had dis- 
cussed transportation to Mars, which would probably 
have been considered ridiculous by the press... .” He 
took certain steps, firmly. On January 19, 1920, in the 
Times, he stated flatly, ““‘Too much attention has been 
concentrated on the proposed flash powder experiment, 
and too little on the exploration of the atmosphere.” 
He tried, vainly, to instill a scientific attitude in news- 
paper readers: 

The point is this: Whatever interesting saan 
there may be of the method that has been proposed, 
other than the purpose for which it was intended, 
no one of them could be undertaken without first 
exploring the atmosphere. Any rocket apparatus 
for great elevations must first be tested at various 
moderate altitudes....Hence...an investigation 
of the atmosphere is the work that lies ahead. 
(Need a parallel be drawn between this and the early 

shots at Cape Canaveral?) 

After disposing, Goddard proposed raising 50 to 100 
thousand dollars by popular subscription ‘“‘to be used 
by the Smithsonian Institution in preparing for, and 
undertaking, a preliminary exploration of the atmos- 
phere.” Naturally, not much if anything was ever heard 
of the subscription again. 

Although not altogether clear, it would seem that 
Goddard was working on a rocket that was scheduled 
for launch. On April 29, 1920, the Tzmes carried a story 
that the National Geographic Society said the test was 
set for July, the object of which was ‘‘the study of the 
aurora and the short wave lengths in the spectrum of 
the Sun which are entirely absorbed by air.” Yet the 
headline said: 


MOON ROCKET TEST 
IS SET FOR JULY 


The rocket would have been propelled by explosion of 
successive charges. Guidance would have been ob- 
tained “through the use of photosensitive cells which, 
through influence of rays of light, [would] keep the 
rocket on its course.”’ He was drawing on his 1908 
notebook entry. 

Five months later, September 20, 1920, a Page 1 
box announced a test within a month. Goddard was 
quoted as saying that he would attempt no great 
altitude “for the simple reason that the funds avail- 
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able... will permit the completion of only a small 
model of about six pounds weight, with a capacity of 
about 60 charges.... The chief reason why the work 
is proceeding slowly is the lack of adequate suppor ah 

A January 28, 1921, story announced the shot for 
“early next Summer.” 

It never came off. 

His complaint of inadequate support is not borne 
out, for he was busy all this time. In 1920 he obtained 
three important patents, one for “Means for Producing 
Electrified Jets of Gas.”’ From 1920 to 1922, on a grant 
from Clark University, he performed various experi- 
ments with liquid oxygen and liquid hydrocarbons, 
such as gasoline, propane, and ether. These were his 
most important experiments—he was removing himself 
from the solid propellant field—and they were financed. 
He arrived at the conclusion that the most “practical 
combination” was liquid oxygen and gasoline. He 
literally never changed his mind. He utilized a proving 
stand and carried out secret experiments near Auburn, 
Massachusetts. 

In these experiments it was shown that a rocket 
chamber and nozzle... could use liquid oxygen 
together with a liquid fuel, and could exert a lifting 
force without danger of explosion and without 
damage to the chamber and nozzle. These rockets 
were held by springs in a testing frame, and the 
liquids were forced into the chamber by the pres- 
sure of a noninflammable gas. 

Goddard had learned his lesson about the wayward 
press and withheld information for a number of years. 
News of “A Method” got to Germany, naturally, and 
in May 1922 he received a letter in quaint English from 
Heidelberg: 


Dear Sir 

Already many years I work at the problem to pass 
over the atmosphere of our Earth by means of a 
rocket. When I was now publishing the result of my 
examinations and calculations, I learned by the 
newspaper, that I am not alone in my inquiries and 
that you, dear Sir, have already done much important 
works at this sphere. In spite of my efforts, I did not 
succeed in getting your books about this object. 
Therefore I beg you, dear Sir, to let them have me. 
At once after coming out of my work I will be honored 
to send it to you, for I think that only by common 
work of the scholars of all nations can be solved this 
great problem. 


A copy of “A Method” was dispatched immediately. 
The German mathematical student brought out at his 
own expense in 1923 in Miinich Die Rakete zu den 
Planeteniumen (The Rocket into Interplanetary Space). 
Its author was Hermann Oberth. 

The year 1923 also found Professor Goddard ready 
for a liquid-fuel test. On November 1, he static-fired a 
small rocket. It was tied down. Results were encourag- 
ing but he wasn’t satisfied. “For one thing, there was 
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the problem of getting the fuels from the tanks im 
the combustion chamber fast enough. He had use 
small pumps....’’ So almost at once Goddard e# 
countered a problem with liquid fuels that still is m 
solved altogether. He was, step by step, entering} 
alone—the realm of systems engineering and techni¢ 
direction. He had touched upon vehicle development 
and he was now into propulsion. 

Yet he must have perceived success, for he went be 
fore the American Association for the. Advancement | 
Science at the end of December 1923 and told 
meeting at Cincinnati that he had ‘only one mo 
step to make before he could prepare a model for flighty 


4d 


some 50 miles above Earth. (We might ask: Was the 


iC 


i] 
: 
“step” the pumps?) The rocket would carry “ver 
delicate apparatus” and would open and close auta 
matically when the rocket reached a certain altitude 
thereby “locking into a chamber a quantity of th 
air.”’ Upon recovery of the rocket, Goddard would 
amine the specimen in his laboratory. He had som 
unusual ideas about the chemosphere: he was certail 
that frozen nitrogen would be found at 50 miles. 

In addition to the air trap, Goddard’s persond 
space physics department would measure pressure 
electrical effects, radiation, and temperature at various 
altitudes. 

The AAAS meeting was further informed that th 
rocket would be fired from Earth at a speed of 6 mi/ see 
“six times faster than the speed ever attained by | 
cannon ball.’’? Such speed ‘‘would enable the rocket ti 
free itself from the attraction of the Earth and keep ob 
travelling for the desired distance.” It would be p ‘ 
pelled by liquid oxygen. The skin would be highl 
polished and marked ‘‘so that the mee | 


[could] keep it under observation long enough to ca 
culate its range and course and to estimate its landini 
place.” 

This is both amusing and admirable, in light 
present-day state-of-the-art. And it is illustrative of | 
strange quirk: a desire for publicity, after all. It is ver! 
doubtful he seriously believed in such a speed. 

On New Year’s Day, 1924, having returned to Wo 
cester, Goddard delivered himself of a familiar plaint 
“Work on the high altitude rocket must be supporte: 
and models supplied for actual trial flight during thi 
coming year, if America is to continue her lead in thi 
branch of scientific research.” 

Perhaps either to encourage such support or t: 
solicit it, Goddard allowed himself to be the subject o 
a Sunday feature in the New York Times, May 25 
1924. Under the large headline, “HOPES TO REACE 
MOON WITH A GIANT ROCKET,” was a lon 
article about him and his rocket. However, now tha 
he knew the methods of the press, he was very cautiou 
to discuss little but his “giant” experimental rocket 
60 inches long and 6 inches in circumference. 

An implied comparison with M. Jules Verne failed t 
elicit any “wild blue yonder” comments. Instead, God 


dard discussed propulsion: “successive explosions” vs 
regular propulsion, supplied by gases released in a 
teady flow.” He favored the latter. His scientific 
lanner is well put: 

I am working with well-known gases, but have found 
new ways of applying their use. The experimental 
rocket also will have considerable mechanism. I need 
not say that it has required years to study the many 
details involved. Almost every step has been an 
experiment within itself. But I have proved ... that 
the rocket will rise practically any desired distance. 
This being true, there is nothing to prevent the con- 
struction of a greater rocket which would go up to 
the Moon. 

This set off the reporter and he again switched to 
Verne and his cannon. He tried to inveigle Goddard 
into discussing a ‘‘space car,’’ but his effort is wasted. 

' Nevertheless, realization that publicity was not 
without its value must have led Goddard to allow 
resident W. W. Atwood, Clark University, to an- 
ounce at commencement exercises in June “‘the solu- 
on of the problem of an explosive for the Goddard 
ocket. A liquid will be used to give propelling force 
without overheating the projectile.” Atwood added 
that a model of the rocket was complete and would be 
sent up later that Summer. It was never fired. 
_ During these early years of work on liquid-fueled 
rockets, the Smithsonian Institution, under Dr. Charles 
G. Abbott, gave Goddard grants; as a matter of fact, 
from 1922 to 1930. The Carnegie Institution, as well, 
under the direction of John C. Merriam, was interested 
am the work and arranged for a substantial financial 
eontribution to Clark. The school, also, added its own 
aid. 

Firmly battened by such a combination, Professor 
Goddard was able to prepare for the signal moment of 
his career, March 19, 1926. 

There was another moment, a very unscientific 
moment, we must not ignore—his marriage on June 
21, 1924, to Esther Christine Kisk of Worcester. A well 
educated woman she was, having been graduated 
from Johns Hopkins. She met Goddard at Clark and 
became his secretary. After marriage, she was present 
at almost every rocket launch and she was the “‘official”’ 
photographer. 

Late in 1925, the New York Times was tracking 
Goddard again. They ran a story about the fated 
young Austrian pilot-astronomer-rocketeer, Max Valier, 
and his plans to reach the Moon by rocket travel. The 
story appeared in the Sunday Magazine Section, 
November 22 and was 90 per cent concerned with 
Valier; however, at the very end Goddard was inter- 
‘viewed. Instead of praise or enthusiasm, he showed 
stiffness and talked very coldly about experiments. He 
probably had on his mind the worries with his liquid- 
fuel static tests. He was ready to try again after almost 
‘two year’s preparation. He felt he had overcome the 
problem of the pumps. In December he tested his 


second rocket. The fuels were forced “into the chamber 
by the pressure of . . . nitrogen.’”’ It worked well. 

There remained but to fly it. 

March 16, 1926, was cold. Goddard had no classes 
that day. There were two inches of snow on the ground 
and the ground was hard and the grass stubbly. God- 
dard, his wife, and his assistant, Henry Sachs, loaded 
the equipment on a trailer and then climbed in a road- 
ster with the isinglass curtains up. They drove quietly 
and quickly to Aunt Effie Ward’s farm near Auburn, 
Massachusetts. She didn’t mind—anything that 
Robert did was all right with her. 

Goddard and Sachs set up the launching stand and 
attached the rocket. Both men posed by their creation, 
bundled against the cold. Goddard wore a cap, over- 
coat, gloves, and buckled overshoes. Sachs simulated 
igniting the rocket. He wore ear-muffs. Esther snapped 
the pictures. The land was bleak and the snow covered 
the base of the launching stand. 

The rocket on the stand was queer looking. The fuel 
tanks were ‘‘slender tubes, placed one behind the 
other.” The combustion chamber, about the size and 
shape of a 4th of July rocket, and exhaust nozzle were 
well ahead—the opposite of today’s design.—and sup- 
ported on “‘spidery arms which also carried the fuel 
lines.”’ Overall length was perhaps 10 feet, about half 
of which was the rocket proper. ‘‘Pressure to force the 
fuels into the combustion chamber was furnished by 
an outside pressure tank and, after launching, by an 
alcohol heater carried on the rocket.’ 

A mistaken idea that this was improvised has man- 
aged to grow. Actually the design was carefully planned, 
and Goddard thought it advantageous in keeping the 
flame away from the tanks and in producing stabiliza- 
tion. For the latter it was of no value. “This is evi- 
dent,’’ Goddard later commented, ‘‘from the fact that 
the direction of the propelling force lay along the axis 
of the rocket, and not in the direction in which it was 
intended the rocket should travel, the condition there- 
fore being the same as that in which the chamber is at 
the rear of the rocket.” 

Henry Sachs had his blow torch going. It was at- 
tached to the end of a long pole. He stood off to one 
side and touched the flame to the igniter at the very 
tip of the rocket motor. Esther turned the crank of her 
camera to record the action. Flame shot from the 
nozzle. The rocket rose 41 feet, turned over, travelled 
184 feet, reached a speed of 60 miles per hour, and 
stayed in the air a total of 2} seconds. 

In the cabbage patch where it had smashed, the 
three of them picked up every smidgen. They had just 
conducted, successfully and secretly, the first liquid-fuel 
rocket flight in the world. But they were Yankees, too, 
and they saved the pieces. 

Not one whisper reached a newspaper. 

Apparently in the next two years, Goddard fired off 
a number of similar rockets; however, the information 
is sketchy and the designs and results are unknown. 


un 
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He was quite busy, at times, with newspaper and 
periodical publicity during this time. On July 4, 1927, 
the New York Times published a story about an 
“ocean rocket” of his that would carry passengers 
across the Atlantic at “terrific speeds,” much faster 
than Byrd or Lindbergh. Goddard played down its 
immediacy and said that there were many details to 
be worked out. (As a matter of fact, it was not until 
1932 before his final thoughts on the idea were made 
known.) Another Times story must have given him 
some small satisfaction, as a scientist. A Sunday 
magazine story carried the huge headline: “A ROCKET 
AUTO OPENS VISTAS OF STAR VOYAGES.” 
More important, especially to Goddard, was the sub- 
head, which was central to his ideas and preachments: 
“But First Rocket Planes Must Test Air Many Miles 
Up and Leap the Atlantic in an Hour and a Half—A 
Principle of Vast Possibilities.”” Goddard was cited in 
the lengthy piece, along with Valier, Opel, Oberth, and 
Esnault-Pelterie. 

In the Summer of 1928, Scientific American called 
on his good services to comment on an article that 
posited a flight from Earth to Mars, propelled by 
smokeless powder, in 7 months. He did not really com- 
mend the piece, but by an odd circumlocution came up 
with this statement: 

According to recent press notices, interplanetary 
transportation must remain impossible until atomic 
energy can be obtained and controlled. This attitude 
is...much like that of the scientists of thirty years 
ago, who declared that an airplane could not operate 
unless the force of gravity could be neutralized. 

If atomic energy were available, it would be a 
very convenient means of propelling an interplane- 
tary rocket. Atomic energy is not, however, neces- 
sary, as an interplanetary flight is possible with 
means even now at our disposal. This is set forth in 
my article in...1919....I can say that, although 
[the] article ...may read like romance, it is never- 
theless thoroughly scientific, and, while not. telling 
the whole story, it gives a good picture of what an 
interplanetary rocket must be like. 

We come, then, to July 17, 1929. 

The Times felt it good enough news next day to run 
a story on Page 2: “METEOR-LIKE ROCKET 
STARTLES WORCESTER.” Goddard called the shot 
“Test of Medium-Sized Rocket Having High Center 
of Gravity and Low Center of Air Resistance’—a 
singularly unimpressive title. The scene was again 
Effie Ward’s farm. But the equipment was far more 
elaborate. A 40-foot tower (also reported as 60-foot) 
stood in a shallow gulley. An 114-foot rocket, fitted on 
two g-inch rails, rested vertically within the frame of 
the tower. The combustion chamber was located at the 
rear of the rocket. It was the best location, Goddard 
was now convinced, because ‘‘no part of the rocket’? was 
in the “high velocity stream of ejected gases, and 
none of the gases” was directed at an angle with the 
rocket axis. Up a slight rise, perhaps 60 to 75 feet 
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away, stood a wooden shack. Inside a small group ‘ 
experimenters were arranged around an electric det 
nator. The rocket, which weighed 32 pounds empty 
held 14 pounds of gasoline and 11 pounds of liqui 
oxygen. It also carried an aneroid barometer, a the 
mometer, and a camera which was focused on bot 
instruments. It was to be operated by a trip lever # 
the ejection of a parachute. Everything was ready. 

Someone in the wooden shack pressed an electr} 
switch. ‘ 

What happened in the next few seconds is not é1 
tirely clear. 

The rocket shot upward with a tremendous roag 
just how high is disputed. Anywhere from 90 to 30 
feet. Near apogee it exploded violently. The noise wé 
heard for two miles around. The townspeople of Wo 
cester were frightened beyond reason. Even as a pars 
chute carried the empty shell of the rocket gently 1 
Earth, the people went into action. At first, the 
thought a huge meteor had exploded. Other witnesse 
insisted that an airplane in flames had flashed acros 
the horizon and blown up. Two police ambulances rac¢é 
to the scene, ‘looking for the supposed victims.” Af 
airplane took off from the local field and submitted #1 
environs to reconnaissance. Newspaper men appearee 

Goddard pleaded that they not publicize the even 
But freedom of the press will out when to its advantagg 
and of course his pleas were ignored. It was too good] 
story. ‘Moony’ Goddard was at it again. Even befo 
they arrived home, Mrs. Goddard recalls, extras we? 
coming out. The professor was forced to issue a stat 
ment. 

The test this afternoon was one of a long series 

experiments with rockets using entirely new pr 

pellants. There was no attempt to reach the Moc 
or anything of such a spectacular nature. The rocke 
is normally noisy, possibly enough to attract con 
siderable attention. The test was thoroughly satisfa 

tory; nothing exploded in the air, and there was 1 

damage except incident to landing. 

He did not really explain the noise—something will 
everyone seems to have been agreed upon. Regardless 
America’s second liquid-fueled rocket had been fire 
publicly and the newspapers were pleased to call it} 
“Moon” rocket. 

And governmental officials were pleased to ax 
Goddard to explain. He called on the Massachuset} 
state fire marshal. Dissatisfied, the marshal withhe} 
permission to conduct further tests. Goddard the 
called on the military authorities at Camp Devens, ‘ 
miles north of Worcester, and asked them for pen 
sion to carry on his field firings. (He had used the cant 
once before.) The military were afraid of fires ar} 
demurred. Determined, the professor appealed to th 
War Department and was supported by his faithf 
Dr. C. G. Abbott of the Smithsonian. On October 2! 
1929, the War Department granted permission to 
an abandoned farm near the artillery range of Ca 
Devens, said permission subject to the understandir 
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Goddard would take precautions against. fires, 
ights were permitted only after a rain or when snow 
on the ground. 

The tower at Auburn was broken down and trans- 
ated to the Camp. Another precaution was a sheet- 
on windshield erected on three sides of the tower. A 
d 50 feet in front of the tower protected the per- 
nnel. The first static test took place on December 3. 
Meanwhile, a series of newspaper and magazine 
icles on Goddard more or less culminated in an 
cle in the Sunday New York Times, October 13, 
29, in which the professor was revealed to be working 
1 a solar or Sun motor. However, the article went on 
to say, ‘‘there is still no method of turning captured 
lar energy efficiently into propulsive energy, for 
ket propulsion requires that some of the mass of the 
ng craft be left behind in space as well as some of 
e energy, this discarded mass providing the ‘action’ 
to which the ‘reaction’ of the rocket is the twin.” 

_ The conclusion was familiar. “Passenger rockets are 
little likely, however, to be the next use of rocket craft. 
Professor Goddard disclaims immediate interest in this 
passenger problem. What he wants is a rocket able to 
pierce most of the Earth’s blanket of air.’ Popular 
Science Monthly, in an editorial in October, had lauded 
Goddard for such an interest: ‘‘While others have 
talked of shooting rockets to the Moon, Professor 
Goddard’s success thus far proves him the most prac- 
tical of the lot.”’ 

_ The Times story was read by Colonel Charles Lind- 
bergh, then riding the top of a fantastic wave of hero 
worship. He was at the time ‘‘particularly interested in 
any research that would provide accurate information 
about the upper layers of the atmosphere.’ On a 
Sunday afternoon in November, the 24th, the Colonel 
slipped into Worcester and visited Goddard’s laboratory 
for two hours. The professor outlined his work in 
progress to the young flier and ran motion pictures of 
his rocket flights. Lindbergh deliberately avoided re- 
porters on his way out of the city, and Goddard refused 
{o comment on the visit. 

But Lindbergh was probably the most important 
guest of Goddard’s life, for out of the visit to the labora- 
tory were to come grants and interest from the Gug- 
genheim family during the 1930’s. However, that is 
another story for another article. 

Shortly after the static tests got underway at Camp 
Devens, the New York Times (indefatigable, when 
seeking after Goddard) ran yet another long story 
about him and his rockets on Sunday, December 15, 
1929. Written by R. L. Duffus, it was a full page long 
and was the most thorough and accurate account yet 
printed. Both the headline and the story emphasized 
that “this Winter” a “great rocket” would be fired to 
xplore outer space. The tone was sober; the informa- 
Hon was wide. About the only drawback are three in- 
sxedibly poor illustrations—and the fact that the 
rocket was not fired. 


Before we leave Goddard and his achievements up to 
1930, we must discuss another aspect of the man which 
we haven’t had time to consider—his inventions. God- 
dard was an inventor and an active one. In Popular 
Science Monthly, October 1929, he wrote an article 
describing and analyzing a rather unusual device: a 
solar motor. Although dated now and even antique, it 
is of some interest. Goddard claimed an efficiency for it 
of at least 50 per cent, “more than twice that of the 
finest steam turbine electric generating plant in opera- 
tion today.” 

In explaining its operation, Goddard pointed to the 
fact that “instead of the usual kind of boiler, a hemi- 
spherical end piece, made of clear fused quartz, is bolted 
to a hollow body. Water is pumped in and mercury 
under high pressure is sprayed into the water... .” 
Quartz would not crack even when subjected to enor- 
mous temperatures. This hemispherical end piece was 
suspended from the apex of three Duralumin struts at- 
tached toa large (about 20 feet) parabolic mirror which 
focused Sun rays on the quartz “‘so that the point of 
greatest heat was where the spray of mercury formed 
a screen in the water.” 

Let Goddard continue in his own way: 
The water itself, being transparent, will not absorb 
the heat, but the mixture, being opaque to light, will 
absorb all the heat instantly. The resulting high 
temperature will convert the metallic mercury into 
mercury vapor at a correspondingly high pressure, 
and the water will be converted into steam at high 
pressure. 

There were at least two advantages to this design: 
(1) the quartz could be shaped so that every ray of 
light would strike it at right angles and (2) the hottest 
point was inside the boiler instead of outside. 

Goddard believed that the commercial possibilities 
were enormous, especially on farms. (This was, recall, 
before the New Deal and Rural Electrification.) A 
100-foot mirror would produce at least 650 horsepower. 
He also considered its use in aerial navigation as well 
space travel. He wrote: 

Recently another possible use...has been men- 
tioned in connection with my experiments in . . . in- 
terplanetary navigation which the Frenchhave named 
“astronautics.”’ 

In 1919 I suggested the possibility of producing a 
rocket so powerful that it would leave the Earth’s 
surface never to return. Since then there has been 
considerable speculation regarding the use of such a 
rocket in... interplanetary space. It is in this con- 
nection that I have been interested in developing a 
very light and efficient solar engine such as now has 
been designed. 

He envisioned a new propulsion. He was entering 
upon an exciting period. The static tests at Devens had 
commenced. The entire Roswell, New Mexico, opera- 
tion lay ahead. He had 15 years of life left him and 
much to do. 
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Technical Notes 


Alignment Chart for Satellite 
| Orbit Calculations 


Leith Holloway* 


In order to track artificial earth satellites during the first 
few days after launch time when accurate machine-derived 
ephemerides are not yet available, it is helpful to be able to 
compute the theoretical values for the rates of motion of the 
nodes and perigee of the orbit and to determine the satellite’s 
mean height from its period or vice versa. First approximation 
formulas for computing these quantities are given in Reference 
(1), and convenient transformations of them are given below. 
However, these formulas are awkward to use because of non- 
integral powers of quantities appearing in them. Therefore an 
alignment chart such as the one published here is a very con- 
venient means for evaluating these formulas. The results ob- 
tained from the chart are probably as accurate as the neces- 
sary approximations justify. See Figure 1. 

Because of the nonspherical shape of the equipotential sur- 
faces of the earth’s gravitational field resulting from the 
earth’s oblateness, the nodes and perigee of a satellite’s orbit 
are continuously changing their positions. These motions are 
akin to the precession of a gyroscope. 

With inclination angles less than 90° the nodes of a satellite 
orbit advance westward with respect to the stars, or in the 
direction of decreasing right ascension. With inclinations 
greater than 90° (retrograde orbits) the nodal motion is east- 
ward. The perigee of an orbit having an inclination angle of 
less than 63.4° moves in the direction of the satellite’s motion, 
whereas above this critical angle the direction of the motion 
is the reverse. The alignment chart: estimates these motions 
by solving the following equations: 


P = $4.49 ath, 
Q = —3.1208 X 105 (cos 7)P-7/8 
@ = 3.1208 X 10° (2 — 2.5 sin? 7) P7738, 


where P is the anomalistic period of the satellite in minutes; 
Ry is the equatorial radius of the earth; a, the semimajor axis 
of the satellite’s elliptical orbit; 7, the inclination angle of the 
orbit; and 2 and & are the rates of motion of the nodes and 
perigee respectively in degrees per day where positive results 
indicate motion in the direction of the earth’s rotation for the 
nodes and in the direction of satellite motion for the perigee. 
The mean height is equal to a minus Rz. 

In the derivation of the constants in these formulas, Ry 
was taken to be 3963.3 statute miles and J, the coefficient in 
the earth’s gravitational potential function, was assumed to 
be exactly 0.0016238, a value recently determined in Refer- 
ence (2). The magnitude of J is related to the degree of flatten- 
ing of the earth at the poles. 

Although the above formulas for 2 and @ are valid only 
for circular orbits, the effect of ellipticity is small with low 
eccentricities, amounting to only two per cent for an eccen- 
tricity of one-tenth. The theoretical value of the correction 
factor is the square of the reciprocal of (1 — e?) where e is the 
eccentricity. 


* U.S. Weather Bureau, Washington D.C. 
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An example will illustrate the use of the alignment cha 
Assume that a satellite has a period of 100 minutes, an - 
clination of 55°, and an eccentricity of 0.12. The left se: 
indicates that this period corresponds to a mean height 
473 miles. Lay a straight edge from this point on the left se: 
through the fifty-five-degree mark on the left side of ‘ 
center scale and extending to the long right scale. Read 
first estimate of the nodal motion on this right scale; this 
3.86 degrees per day or 19.4 minutes of local time earlier ea 
day as read from the left side of this same scale. This tit 
scale is merely the negative of the angular motion of the noc 
converted to minutes plus 3.94, the daily motion of the s 
in minutes. The small correction for the eccentricity of 0.12 
determined from the short scale on the far right of the cha 
The scalar distance given by this short scale corresponds 
the magnitude of the correction to the nodal motion given | 
the long right scale. This distance should be laid off on t 
long right scale upward from the first estimate to the e 
rected value of 3.97 degrees per day. The same procedure 
used to obtain the perigee motion except that the fifty-fiy 
degree mark on the right side of the center scale is employ 
instead. This gives values of 2.17 and 2.23 degrees per d. 
before and after correction respectively. 

It is impractical to extend the scale for the inclinati 
angle used in computing the perigee motion above 60° becat 
of the reversal in direction of this motion at an inclination 
63.4°. Rates of motion of the perigee of high-inclination orb 
must be computed by hand using the formula. In the spee 
case of a polar orbit (¢ = 90°) the magnitude of the perig 
motion is one fourth of what its value would be for an equ 
torial satellite with the same mean height and orbit ecce 
tricity. 

There is an obvious connection between the time scale 
the chart and the prime sweep interval which is the time 1 
quired for the orbit to make one apparent revolution arow 
the earth. Since the orbit is generally precessing towards t 
west against the rotation of the earth, the prime sweep inter\ 
is most often less than one day. In fact, it is usually express 
as one day minus a certain number of minutes. One day min 
the value obtained from the alignment chart’s time scale 
approximately equal to the prime sweep interval. Howev: 
the number of minutes to be subtracted will always be t 
large by about one per cent at 15 minutes and by about t' 
per cent at 30 minutes. The reason for this discrepancy 
that the chart computes the daily rate of the orbit’s local ti 
change, whereas the subtrahend of the prime sweep inter 
expresses this change in units of local time change per peri 
equal to the prime sweep interval itself. 

Clearly the procedure of specifying the position of the orl 
in respect to a coordinate system rotating with the mean s 
has obvious advantages. Knowing the local times under ¢ 
orbit at a given latitude greatly simplifies the problem of fo 
casting the time when the orbit will enter or leave the eart 
shadow at this particular latitude, and it also facilitates p 
dicting when the orbit will be favorably located with resp¢ 
to the twilight zones. This is the justification for including t 
time scale on the chart opposite the angular velocity scale 
nodal motion. Incidentally, this time scale, of course, has 
meaning with reference to the perigee motion and go shot 
be disregarded in this case. 

The author is indebted to all of the authors of Referet 
(2) below for their kind assistance in the preparation of t 
chart and to Mr. Earl Rayfield for assisting with the drafti 
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leferences Orbit of Satellite 1958 6.’’, Journal of Geophysical Re- 
search, Vol. 64, No. 2, February 1959, pp. 209-16. 
1) Myron Lecar, JoHN SoreENsoN, AND ANN Ecxkexs, “A (2) Joun A. O’Keern, ANN Eckets, AND R. KENNETH SQUIRES, 
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Radiation Shelters For Lunar 
Exploration 


Donald H. Robey 


Abstract 


The hazards to human life from solar flare and cosmic 
ray particles have been briefly considered. It is con- 
cluded that human explorers should prepare emergency 
flare shelters with thicknesses of the order of 50 to 100 
or more centimeters of Moon material. This is assuming 
a specific gravity of roughly two for the surface material. 
Reasons are also presented for suggesting that the dark 
side of the Moon may not be entirely safe from flare 
particle radiation. 


Introduction 


Although the actual, manned lunar station is several years 
away, plans are being made for this gigantic undertaking at 
the present time. It is therefore not too early to speculate on 
the nature of shelters required for station personnel. The 
period of time that they remain at the base depends upon the 
radiation environment, the shielding, and the allowable radia- 
tion dose. 

There are two well known active sources of radiation: (1) 
the galaxy and (2) the Sun. The Earth may even be a minor 
source of particle radiation. For example, secondary neutrons 
from cosmic primaries may leave the atmosphere in all 
directions. 

The galaxy is believed to be the source of cosmic radiation. 
About 80 per cent of this radiation is believed to consist of 
protons, 19 per cent of alpha particles (helium nuclei) and the 
remaining | per cent of heavier nuclei such as carbon, nitro- 
gen, oxygen and iron. Electrons are also believed to be present. 

Since cosmic rays are charged, they cannot pass undis- 
turbed through the Earth’s magnetic field unless they are 
traveling parallel to the flux lines. The north and south geo- 
magnetic poles have 90° dips at 76° N, 102° W and 68° §, 
145° K, respectively. Particles incident vertically near these 
poles are not deflected significantly. Fortunately, the atmos- 
phere acts as a radiation shield. In order for a proton to reach 
sea level, it must have a kinetic energy of about 2 Bev in the 
polar regions. However, over the geomagnetic equator, be- 
cause of the EKarth’s magnetic field, the minimum required 
kinetic energy is about 15 Bev. Thus, the equatorial magnetic 
field is equivalent to about 6 atmospheres. Protons in the 
60 Bey range, however, can penetrate the atmosphere, to sea 
level, at all angles of incidence. 

On the surfaces of the Earth, the atmosphere provides us 
with about 1023 g/cm? of shielding against incident particle 
radiation. If the bulk density of the Moon’s surface is taken 
to be about 2 g/cm?, then roughly 511 cm, or more (since the 
range of protons in maria, or rock, for example will be greater 
than in air) of moon material will be needed in place of the 
Harth’s atmosphere. The ratio of ranges for proton penetra- 
tion in almost all of the elements (excluding hydrogen) is less 
than two, so that a more conservative estimate might require 
up to 1023 em (= 34 ft) of moon for the depth of a lunar 
surface shield. 

The additional protection offered by the Earth’s magnetic 
field except for over the polar regions has not been included. 
This could imply an additional 30 or 40 meters of moon. 

The cosmic ray background is believed to be fairly steady 
except for the Forbush intensity decreases which are associ- 
ated with transient activity on the Sun, i.e., solar flares. Since 
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a decrease in cosmic ray background is nothing to be alar ned 
at, it is probably feasible to estimate the maximum radiat: 109 
dose from this source in advance. : | 
A fairly serious hazard from cosmic rays 18 the possible int 
jurious effects on the brain and eyes, especially from heay, 
primaries. According to Davis (1960), some radiobiologis 
believe that single, heavy cosmic particle could incapacitate 4) 
man if it struck in a critical brain area. He estimates that th) 
particles might damage a region of tissue with a diamete? 
equal to about 20 blood corpuscles. If this turns out to be th 
case it is a serious blow to the manned exploration of th 
Moon’s surface, and to space flight in general. It is proba bly. 
too early to properly evaluate this potential hazard at the 
present time. Also the brain can repair itself or use othe? 
channels for many kinds of injuries. ; 
If we take an optimistic attitude and assume that a m at 
could be exposed to this radiation for several days at a ti ej 
and have a good chance of not suffering serious damage, he 
still must cope with solar flare particles. It is known fro nt 
satellite and balloon studies that at least some of the larges 
class 3+ solar flares eject beams of protons, helium nuclei, ane 
electrons, in quantities that are lethal to unshielded man. The 
intensities may be as high as 10* or 104 times the normal cosmi i 
ray background for periods of several hours. The energies are 
ordinarily below the cosmic ray range, i.e., below a few 
hundred Mev, but the flux densities can be orders of magmi 
tude greater. 
However, on eight occasions during the last two decades the 
Sun has erupted and thrown out a burst of high-energy parti! 
cles that might properly be called ‘“solar-cosmic rays.” The 
largest of these occurred on February 23, 1956 and accelerate¢ 
an estimated 10 particles from thermal energies to energie? 
of over 2 Bev (Parker 1957). i 
Since the time of this last great flare, approximately 25 ad. 
ditional events of a similar, but less intense, nature have oc# 
curred. The particles are believed to consist largely of pro‘ 
tons. Furthermore, since the energies and flux densities 0 
these particles can be dangerously high for many hours follow’ 
ing flares, it will be necessary to introduce emergency shield) 
from Moon explorers. | 
When on a field trip several hours distance from a base, 1) 
might be necessary to dig deep holes in the Moon for flare 
shelters. Boards could be carried along to support vehiclet 
and other available equipment over openings to shelters. 
It was estimated (Robey, 1960) that the large flare 0} 
May 10, 1959 would have delivered a radiation dose of abou; 
8000 rem (based on a relative biological effectiveness of twé 
for protons below 40 Mev) to an unshielded man in space aj 
the Earth’s distances. The dose to a man on the Moon’s sur} 
face is about half of this value, if the majority of the particles 
are omnidirectional. 
Actually there is no evidence to suggest that flare particles 
travel in straight paths to the Earth. In the first place, they 
are charged and may be forced to depart from the Sun ix} 
corkscrew orbits depending on the shape of the solar magnetid 
field. In addition to the possible solar dipole field there ard 
reasons to believe (Alfven, 1954) that magnetically polarized 
plasmas are ejected from the Sun. Events of this kind tend té 
create magnetic flux tubes which curve out away from tha 
Sun, because of its rotation, and clear the way for subsequent 
charged particles. Alfven (1954) postulated that fluxes ad 
high as 10-4 gauss may exist in these regions at the Earth’s 
distance. This was subsequently verified by measurements 
made by Pioneer V. Charged particles which follow the fluy 
lines in these tubes might tend to be confined, which coulc 
explain why there does not seem to be a gradual buildup of the 
particles arriving at the Earth. That is, there seems to be é 
sharp leading edge. | 
On the other hand the decay of the particle influx is rela 


vely slow, sometimes lasting for several days. The lifetime 
f the flare, as indicated by the intense emission of hydrogen 
light, ordinarily does not exceed a few hours. If this also 
esponds to the duration of the period of intense particle 
tion, then the particles must be stored somewhere. 
ckler (1959) states that a storage in the vicinity of the 
fits the experimental facts very well, that is, the flare 
icles observed terrestrially, decay slowly, and stream in 
the poles with high intensities both prior to and following 
2 associated Forbush decrease in cosmic ray intensity. The 
dden shock-like arrival of particles might be explained by 
he tendency to shove and carry along quantities of ionized 
terplanetary ions. This is analogous to the action of the solar 
nd on ionized nitrogen in comet tails. The tails are acceler- 
ed and blown out from the Sun with forces greatly exceeding 
he radiation pressure of sunlight. 

Meyer, Parker and Simpson (1956) argue in favor a field- 
e (< 10~® gauss) heliocentric cavity extending out to about 
astronomical units (AU’s). The lack of significant mag- 
tic field in this region can be used to explain the sharp onset 
osmic noise absorption over the poles (attributable to the 
rival of protons, which increases the ionization in the 
)-layer and, consequently, attenuates radio noise from space). 
Phere is evidence to suggest that the flux is omnidirectional. 
th this circumstance as well as the observed slow decays in 
ensity can be explained by a scattering and diffusing prop- 
arty of the boundary of the heliocentric cavity. The boundary, 
ghich is actually apt to be a spheroidal shell, is probably 
yeral AU’s thick. That is, it is probably in a turbulent or 
li sordered state in a shell of thickness comparable to the or- 
yital distances of Mars and Jupiter. The Sun is believed to 
teate this cavity by the continuous action of the solar 
vind (Parker, 1958). Also the occasional ejections of mag- 
etically polarized plasmas are probably able to modulate 
he magnetic properties in the boundary layer (between 
Mars and Jupiter) and thereby bring on Forbush decreases. 
_ The implication of a heliocentric cavity is that the flare 
articles can scatter back and forth in the solar system and 
roduce an omnidirectional flux. Consequently, a colony on 
he dark side of the Moon would not be safe from this radia- 
ion. 

If the particles leave the Sun in a horn shaped beam, as 
uggested in Fig. 1, and are scattered at the interstellar mag- 
etic field boundary, it is reasonable to suppose that there 
vill be a tendency for them to remain trapped within this 
listurbed region created by the original beam. This prevents 
hem from scattering immediately all over the inner solar sys- 
em. Particles with higher energies are expected to diffuse 
ut of the beam at faster rates. 

If the particles were not restricted by some mechanism 
uch as this, flares on the back side of the Sun would produce 
ffects measurable at the Earth, e.g., cosmic noise absorption 
ver the magnetic poles. However, the flux rate of growth 
vould be gradual and probably less intense by an order of 
pagnitude, rather than a sharp spike which builds up in 
ainutes. We would expect the number of gradual buildups to 
i comparable to the number of sharp increases. However, 
he gradual buildups apparently are not observed except 
or rare cases, where large flares occur within two or three 
ays after the region passes over the west limb. 


we 


ummary 


Since the Moon lacks a natural radiation shield such as 
n atmosphere, or a significant magnetic field, it will be neces- 
ary for human explorers to have radiation shelters. 

It is suggested that the living quarters and working areas 
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Fic. 1 showing the tendency for solar plasmas to follow 
magnetic flux lines. The high energy proton beam from the 
flare of 2/23/1956 may have missed the Earth. The first pro- 
tons to arrive at the Karth apparently leaked out of the beam 
as indicated by arrows. 


be located beneath the surface, or inside a mountain, so that 
at least 1023 gm/cm? of material is between the personnel and 
the incoming radiation. This roughly duplicates the shielding 
from the Earth’s atmosphere but neglects the additional pro- 
tection from the Earth’s magnetic field. However, there is 
only about a 10 per cent difference in cosmic ray background 
at sea level in going from the geomagnetic equator to a geo- 
magnetic latitude of about 40°. The changes are small there- 
after. 

When making manned explorations on the surface it is ad- 
visable to have emergency shelters available because of the 
hazard from solar flares. Directed, rocket-borne explosives, 
e.g., shaped charges, might be used to create hole shaped 
shelters. About 50 or 100 em of Moon material would proba- 
bly shield most flares. The personnel might have to remain in 
the shelter for several hours. The rarer, cosmic ray-type flares, 
may require two or three meters of Moon. 

Reasons are given for expecting an approximately isotropic 
distribution of the particle influx on the Moon. It is concluded 
that explorers on the dark side will therefore also be in danger 
from solar flare radiation. 

A very small, but extensive, lunar magnetic field might 
offer some protection. 
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Greek Alphabet 


A a alpha = (a) N yp nu (n) 
B £6 beta (b) = &€ xi (x) 
T y gamma (g) O o  omicron (6) 
A 6 delta (d) Iz pi (p) 
E ¢« epsilon (@) P p_ rho (r) 
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6 6 theta (th) YT v_ upsilon (u) 
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The Orbital Elements 
a = semi-major axis 
e = eccentricity 
§3 = longitude of the ascending node 
i = inclination to plane of the ecliptic 
w = longitude of perihelion measured from the 
node 
T = time of perihelion passage 
Complicated exponents and subscripts should be 
avoided when possible to represent by a special symbol. 
Fractions in the body of the text and fractions occur- 


ring in the numerators or denominators of fractions 
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cos (7x/2b) 
cos (7ra/2b) 


is the preferred usage. 

The intended grouping of handwritten formulas can 
be made clear by slight variations in spacing, but this 
procedure is not acceptable in printed formulas. To 
avoid misunderstanding, the order of symbols should 
therefore be carefully considered. Thus: 

(a + bx) cos t is preferable to cos t (a + bx) 

In handwritten formulas the size of parentheses, 
brackets and braces can vary more widely than in 
print. Particular attention should therefore be paid to 
the proper use of braces, brackets, and parentheses 
(which should be used in this order). Thus: 


{[a + (b + ex)*] cos ky}? 
is required rather than ((a + (b + cx)") cos ky)?. 


Equations are numbered and referred to in text as 
(15). 
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